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geometric  variables. 

The  nozzle  and  boattail  contours  are  modeled  as  a second-order  polynomial 
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—and  a cone,  respectively.  For  fixed  initial  expansion  contours  and  nozzle 

and  boattail  lengths,  the  second-order  polynomial  nozzle  contour  is  uniquely 
specified  by  the  nozzle  throat  attachment  angle  and  the  nozzle  exit  lip 
radius,  and  the  conical  boattail  is  uniquely  specified  by  the  boattail  exit 
lip  radius.  These  three  independent  parameters  are  varied  to  determine  the 
unique  nozzle-base-boattai 1 configuration  that  yields  maximum  thrust. 

Three  methods  are  included  in  the  program  to  determine  the  maximum  thrust 
contour.  Each  method  requires  an  initial  estimate  of  the  geometry  that  pro- 
duces maximum  thrust.  The  methods  are  derivative  methods  that  perturb  the 
wall  geometry  paranietVs  to  determine  approximate  local  derivatives.  This 
informaticn  is  used  tori  calculate  a search  direction  and  a step  length  in  the 
optimization  variable  space.  Successive  applications  of  the  procedure  modify 
the  contour  parameters' until  the  maximum  thrust  configuration  is  obtained. 

The  first  method  is  axial  iteration.  It  treats  the  thrust  as  a function 
of  up  to  three  independent  coordinates,  and  moves  in  the  direction  of  each 
coordinate  individually,  using  first  and  second  partial  derivative  approxi- 
mations to  determine  the  step  length.  The  second  method  is  a gradient  tech- 
nique, the  method  of  steepest  descent.  Jt  computes  an  approximation  to  the 
gradient  of  the  thrust  function  and  steps  an  arbitrary  step  length  in  the 
direction  of  the  gradient.  The  third  technique  is  Newton's  method,  a second- 
derivative  method.  Newton's  method  computes  first  and  second -derivatives  of 
the  thrust  function,  defining  both  the  direction  and  the  step  length  for  the 
optimization  step. 

The  three  methods  were  first  applied  to  the  design  of  maximum  thrust 
nozzle  contours.  Each  of  the  methods  converges  efficiently  to  the  maximum 
thrust  nozzle  contour.  Newton's  method  and  axial  iteration  converge  very 
quickly,  in  approximately  twenty  nozzle  flow  field  calculations  or  less. 

The  method  of  steepest  descent  requires  as  much  as  twice  as  many  flow  field 
calculations  as  the  other  methods,  but  proceeds  to  the  exact  maximum  by  con- 
verging to  a point  whose  gradient  is  zero.  The  nozzle  contours  from  these 
methods  are  shown  to  yield  nearly  identical  thrusts  to  those  predicted  by 
maximization  using  the  calculus  of  variations,  affirming  the  validity  of 
the  direct  optimization  approach. 

”^The  three  methods  were  also  applied  to  the  design  of  nozzle-base-boattail 
configurations.  The  three-dimensional  optimizations  were  found  to  be  prac- 
tical for  finding  the  maximum  thrust  attainable  for  a nozzle-base-boattail 
assembly.  Each  of  the  three  methods  converged  efficiently  to  the  maximum 
thrust  nozzl e-base-boattail  contour. 
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SECTION  I 


INTRODUCTION 


Current  design  techniques  for  propulsive  nozzle  contours  assume 
that  the  nozzle  exhausts  into  a static  region  having  a specified 
ambient  pressure.  The  design  objective  is  to  select  that  contour  which 
yields  maximum  thrust  for  that  ambient  pressure.  As  illustrated  in 
Fig.  1,  the  thrust  producing  region  is  much  more  complex.  The  nozzle 
is  generally  installed  in  an  airframe  with  an  annular  base  region  be- 
tween the  nozzle  lip  and  the  airframe  lip.  Often  the  aft  portion  of 
the  airframe,  referred  to  as  the  boattail , is  contoured  to  lessen  the 
thrust  loss  associated  with  the  low  pressure  base  region.  Considera- 
tion of  boattail  and  base  thrusts  modifies  the  optimum  nozzle  contour 
from  that  for  an  isolated  nozzle. 

A technique  for  designing  conical  nozzles  in  a cylindrical  boat- 
tail  including  the  effects  of  the  base  region  was  developed  by  Byington 
and  Hoffman  (1),  making  possible  performance  increases  of  up  to  two 
percent.  No  analogous  technique  is  available  for  designing  contoured 
nozzles  including  the  effects  of  the  base  region  or  for  designing  con- 
toured nozzles  and  contoured  boattail s including  the  effects  of  the 
base  region. 

The  objective  of  the  present  investigation  is  to  develop  an  effi- 
cient method  for  the  design  of  maximum  thrust  contours  for  nozzle- 
base-boattail  assemblies,  such  as  illustrated  in  Fig.  1. 

The  nozzle  and  boattail  contours  must  satisfy  certain  geometric 
and  physical  constraints;  for  example,  fixed  overall  length  for  both 
contours,  specified  throat  radius  and  initial  expansion  contour  for 
the  nozzle,  specified  outer  radius  and  initial  expansion  contour  for 
the  boattail,  and  specified  initial  conditions  for  the  nozzle  and 
boattail  flow  fields.  Substantial  effort  has  been  directed  to  proce- 
dures for  maximum  thrust  nozzle  design,  and  several  computer  programs 
for  nozzle  design  based  on  the  calculus  of  variations  have  been  de- 
veloped. The  difficulty  with  the  calculus  of  variations  approach  is 
that  if  the  gas  dynamic  model  is  modified  or  the  general  geometry  of 
the  thrust  producing  region  is  changed,  the  entire  optimization  anal- 
ysis and  computer  program  must  be  reworked,  a task  requiring  several 
man  years  of  effort. 


Byington,  Jr.  and  J.D.  Hoffman,  "Effects  of  Base  Pressure  on 
Conical  Thrust  Nozzle  Optimization,"  AIAA  Journal,  Vol . 7,  No.  3, 
March  1970,  pp.  380-382. 
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Figure  1.  Combined  Nozzle-Base-Boattail  Optimization  Model 
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Tlte  approach  taken  in  Uiis  investigation  is  to  develop  an  opti- 
mization procedure  that  is  valid  <^or  any  gas  dynamic  model  and  a 
variety  of  geometric  configurations.  To  achieve  that  goal,  the  nozzle 
wall  contour  was  assumed  to  oo  a second-order  polynomial  of  the  form 
y = a + bx  + cx^,  and  the  boattail  contour  was  assumed  to  be  a first- 
order  polynomial  of  the  form  y = a + bx.  Nonlinear  programiiing  methods 
(2)  are  employed  to  vary  the  polynomial  coefficients  a,  b,  and  c, 
seeking  the  nozzle  and  boattail  contours  that  yield  maximum  thrust. 

Two  questions  must  be  answered.  First,  do  wall  contours  specified 
by  polynomials  yield  results  comparable  to  those  obtained  by  the  calcu- 
lus of  variations?  Second,  which  nonlinear  programming  methods  most 
efficiently  solve  the  nonlinear  programming  problem? 

The  present  analysis  is  limited  to  applications  where  the  external 
flow  over  the  boattail  is  supersonic.  Viscous  effects  are  neglected 
in  both  the  nozzle  flow  field  and  the  boattail  flow  field.  Consequently, 
E transonic  flow  effects  and  flow  separation  effects  on  the  boattail 

1 are  both  neglected.  The  resulting  computer  program  is  thus  limited  to 

supersonic  flight  speeds  at  design  conditions  where  separation  does 
i not  occur.  The  optimization  procedure,  however,  does  not  depend  on 

[ the  flow  field  model,  and  it  should  be  possible  to  apply  the  procedure 

j to  any  flight  speed,  including  viscous  effects,  when  rapid  and  accurate 

numerical  methods  for  determining  the  cor>^e5ponding  flow  field  become 
available. 


I 

r 

i 

i 

i 


i 

I 

i 


— 

' 'D.M.  Himmelblau,  Applied  Nonlinear  Programming,  McGraw-Hill  Book 
Co.,  New  York.  1972. 
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SECTION  II 


ANALYSIS 


The  design  of  maximunt  thrust  contours  for  any  thrust  producing 
element  requires  the  ability  to  analyze  the  flow  field  in  that  partic- 
ular element,  and  an  optimization  method  for  determining  the  maximum 
thrust  contour.  The  flow  field  analysis  procedure  and  the  optimization 
methods  employed  in  this  investigation  are  discussed  in  this  section. 


1 . Flow  Field  Analysis 

The  nozzle-base-boattai 1 assembly  geometry  is  illustrated  in  Fig.  1. 
The  initial  expansion  contours  for  both  the  nozzle  and  the  boattail  are 
assumed  to  be  circular  arcs.  The  nozzle  contour  is  assumed  to  be  a 
second-order  polynomial  of  the  form 

y(x)  = a + bx  + cx^  (1) 

and  the  boattail  contour  is  assumed  to  be  a first-order  polynomial  of 
the  form 


y(x)  = a + bx  (2) 

The  above  choices  for  the  nozzle  and  boattail  contours  require  the 
determination  of  the  coefficients  a,  b,  and  c for  the  nozzle  contour 
and  the  coefficients  a and  b for  the  boattail  contour.  These  coeffi- 
cients are  uniquely  determined  by  specifying  the  nozzle  throat  attach- 
ment angle  ©an.  the  nozzle  exit  radius  yen>  the  boattail  exit  radius 
yeb»  and  by  requiring  that  the  polynomial  contours  attach  smoothly  to 
the  circular  arc  initial-expansion  contours  (see  Fig.  1).  These  three 
parameters  uniquely  specify  a nozzle-base-boattail  configuration.  The 
optimization  procedure  searches  through  the  allowable  ranges  of  these 
three  parameters  to  determine  the  particular  set  of  values  that  yields 
the  maximum  total  thrust  on  the  nozzle-base-boattail  assembly. 

The  gas  dynamic  model  is  that  for  the  steady  axi symmetric  flow  of 
an  inviscid  fluid  in  the  absence  of  work,  heat  transfer,  and  body 
forces.  Such  a flow  is  isentropic.  The  governing  partial  differential 
equations  are  the  following  (3). 


pu^  + pVy  + up^  + vpy  + 6pv/y  = 0 

(3) 

puu^  + pvu^  Px  ^ ° 

(4) 

4 


puv^  + pvv^  + p^  = 0 

(5) 

2 2 

up  + vp  , a Up  - a vp  = 0 
^ X y X ' y 

(6) 

Equation  (3)  is  the  continuity  equation,  equations  (4)  and  (5)  are  the 
X and  y Euler  momentum  equations,  and  equation  (6)  is  a form  of  the 
speed  of  sound  equation  for  an  isentropic  flow. 

This  system  of  four  partial  differential  equations  may  be  replaced 
by  a system  of  characteristic  and  compatibility  equations,  which  are 
ordinary  differential  equations  (see  Appendix  A).  For  isentropic  flow, 
the  acoustic  speed  is  a function  of  the  local  static  pressure  and 
density.  The  composition  is  assumed  to  be  a general  mixture,  either 
fixed  (frozen)  in  composition,  or  in  chemical  equilibrium.  The  gas 
dynamic  model  requires  that  the  entropy,  the  stagnation  pressure,  and 
the  stagnation  enthalpy  remain  constant  along  streamlines,  although 
these  properties  may  vary  between  streamlines.  A conventional  predictor- 
corrector  second-order  numerical  integration  procedure  (4)  is  employed 
to  integrate  the  characteristic  and  compatibility  equations,  starting 
from  a supersonic  initial-val ue  line. 

Several  procedures  have  been  developed  for  analyzing  the  transonic 
flow  in  the  throat  region  of  convergent- di vergent  nozzles.  Kliegel's 
method  (5)  is  employed  in  the  present  study  because  it  is  accurate  and 
easy  to  use.  Kliegel's  method  stems  from  a study  by  Hall  (6).  Both 
solutions  involve  expansions  for  the  velocity  components  u and  v in 
inverse  powers  of  R,  the  ratio  of  the  upstream  throat  wall  radius  of 
curvature  ptu  to  the  throat  radius  yt  (R  = Ptu/yt)-  Rail  used  a 
straight  expansion  in  inverse  powers  of  R.  However,  the  resulting 
series  is  divergent  for  R < 1.  Kliegel  modified  Hall's  approach, 
developing  an  expansion  in  inverse  powers  of  (R  + 1),  which  is  conver- 
gent for  R < 1.  The  complete  third-order  axisymmetric  solution  for 
the  velocity  components  derived  by  Hall  and  Kliegel  is  presented  in 
Appendix  B.  Thus,  the  velocity  components  along  an  initial-value  line 
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may  be  calculated,  defining  a starting  line  for  the  numerical  solution 
of  the  supersonic  portion  of  the  nozzle  flow  field  by  the  method  of 
characteristics . 

The  nozzle  flow  field  is  illustrated  in  Fig.  2.  Contour  AA'  is 
the  circular  arc  initial  expansion  contour.  Point  D is  the  nozzle 
throat  attachment  point  where  the  maximum  thrust  contour  DF  attaches 
snoothly  to  the  initial  expansion  contour.  The  maximum  thrust  nozzle 
contour  is  obtained  by  varying  both  the  location  of  point  D along  the 
initial  expansion  contour  and  the  contour  between  points  D and  F.  The 
initial-value  line  AB  is  obtained  by  applying  Kliegel's  method,  or  by 
defining  the  line  point  by  point  from  any  transonic  flow  analysis. 
Starting  from  the  initial-value  line,  the  method  of  characteristics  is 
applied  to  calculate  the  flow  field  in  region  ABCFA.  Region  ABCEDA  is 
termed  the  kernel,  and  is  affected  solely  by  the  initial-value  line 
properties  and  the  nozzle  i ni tial -expansion  contour.  Region  DEFD,  which 
is  influenced  by  the  nozzle  wall  contour,  is  the  flow  field  region 
bounded  by  the  right-running  Mach  line  DE  originating  on  the  wall  con- 
tour at  the  attachment  point  D and  the  left-running  Mach  line  EF  passing 
through  the  nozzle  exit  lip  point  F. 

The  analysis  program  employed  in  the  study  is  a modification  of 
the  program  developed  by  Pasley  and  Hoffman  (7).  That  program  was 
modified  to  minimize  the  computational  effort  during  optimization,  where 
many  similar  flow  fields  are  computed.  This  is  particularly  pertinent 
in  the  nozzle,  where  much  computational  effort  is  spent  calculating  the 
flow  field  in  the  kernel.  In  general,  the  nozzle  throat  attachment 
angle  0an  corresponding  to  point  D has  some  minimum  value  below  which 
the  optimization  search  will  never  proceed.  The  portion  of  the  kernel 
upstream  of  that  point  is  fixed  and  need  not  be  recalculated  as  point  D 
and  the  contour  DF  are  varied  during  the  optimization  search.  That 
minimum  nozzle  throat  attachment  point  is  indicated  as  point  G on  Fig.  2. 
Consequently,  right-running  Mach  line  GH  may  be  stored  and  employed  as 
the  initial-value  line  for  subsequent  nozzle  flow  field  calculations. 

The  secondary-start  line,  as  this  new  initial-value  line  is  called,  may 
save  as  much  as  50  percent  of  the  computational  effort  required  for  a 
nozzle  flow  field  analysis.  Care  must  be  exercised,  however,  not  to 
compute  the  secondary  start  line  at  too  large  an  angle,  since  subsequent 
nozzle  flow  field  calculations  must  have  throat  attachment  angles 
greater  than  that  of  the  secondary-start  line.  The  program  is  not 
capable  of  restarting  at  any  nozzle  throat  attachment  angle  less  than 
or  equal  to  that  of  the  secondary-start  line. 

The  left-running  Mach  line  Id  is  the  initial-value  line  for  the 
boattail  flow  field  analysis.  All  of  the  flow  properties  along  the 
ini tial -val ue  line  must  be  known  from  a prior  analysis  of  the  flow 
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around  the  fonvard  section  of  the  vehicle.  If  such  information  is 
unavailable,  a uniform  flow  i ni t la 1 -val ue  line  may  be  employed. 

Contour  II'  is  the  circular  arc  initial  expansion  contour,  and  point  L 
is  the  point  where  the  conical  boattail  contour  KL  attaches  to  the 
initial  expansion  contour  II'.  The  external  flow  field  is  assumed  to 
be  supersonic,  so  that  the  iDethod  of  characteristics  may  be  employed 
to  calculate  the  flow  field  in  region  IJKI.  That  region  is  bounded 
by  the  left-running  Macli  line  Id,  which  is  the  initial-value  line,  and 
the  right-running  Mach  line  JK  passing  through  the  boattail  exit  lip 
point  K. 

The  boattail  contour  and  nozzle  contour  may  have  a base  region 
between  them,  as  illustrated  in  Fig.  2.  The  total  thrust  on  the  vehicle 
must  account  for  the  base  pressure  in  this  annular  region.  Several 
models  have  been  developed  for  predicting  the  base  pressure  aft  of  a 
supersonic  nozzl e-boattai 1 assembly.  The  empirical  base  pressure  model 
developed  by  Addy  (8)  was  employed  in  the  present  investigation.  A 
curve  fit  of  Addy's  data,  similar  to  that  employed  by  Byington  and 
Hoffman  (1),  was  employed  to  relate  the  base  pressure  Pb  to  the  exit 
conditions  of  the  nozzle  and  boattail.  This  curve  fit  has  the  follow- 
ing form. 

The  four  r,  functions  are  correlation  factors,  which  are  applied  to  the 
nominal  case  for  which  f(Pen/Peb)  was  obtained.  Each  factor  accounts 
for  variations  in  the  parameter  denoted  within  the  corresponding 
parentheses.  The  nominal  case  and  the  correlation  factors  are  discussed 
in  Appendix  C.  Correlation  factors  based  on  specific  heat  ratio  y and 
the  base  region  radii  yen  and  yeb  have  been  neglected,  since  their 
influence  on  the  base  pressure  is  slight. 

The  total  thrust  of  the  nozzle-base-boattai 1 assembly  is  the  sum 
of  the  thrusts  of  the  nozzle,  the  base,  and  the  boattail. 

The  program  developed  during  this  investigation  can  perform  the 
following  types  of  calculations. 

1.  Analysis  of  a specified  nozzle-base-boattail  assembly. 

2.  Design  of  a maximum  thrust  nozzle  contour. 

3.  Design  of  a maximum  thrust  boattail  contour. 

4.  Design  of  a maximum  thrust  nozzle-base-boattail  contour. 


^^^A.L.  Addy,  "Analysis  of  the  Axisymmetric  Base-Pressure  and  Base- 
Temperature  Problem  with  Supersonic  Interacting  Freestream-Nozzle 
Flows  Based  on  the  Flow  Model  of  Korst,  et  al..  Part  I:  A Com- 

puter Program  and  Representati ve  Results  for  Cylindrical  After- 
bodies," Advanced  Systems  Laboratory,  Redstone  Arsenal,  Alabama, 
Report  No.  RD-TR-69-12,  July  1969. 
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2.  Optimization  Procedures_ 

The  nozzle  and  boattai I flow  fields  are  entirely  independent,  and 
are  calculated  separately.  Thus,  the  nozzle  and  boattail  contour 
parameters  are  related  only  when  considering  a nozzle-base-boattail 
assembly,  because  they  influence  base  pressure  and  hence  total  thrust 
on  the  assembly. 

The  nozzle  typically  provides  the  greatest  thrust  contribution  to 
the  vehicle.  The  boattail  and  base  regions  are  normally  contoured  in 
such  a way  as  to  minimize  the  low  pressure  base  region  and  boattail 
divergence  loss  effects. 

The  approach  taken  in  this  investigation  is  to  model  the  wall 
contour  geometry  as  a second-order  polynomial  for  the  nozzle  contour 
and  a first-order  polynomial  for  the  boattail  contour  (three  independent 
parameters).  These  parameters  are  systematically  varied  to  find  the 
maximum  thrust  nozzle-base-boattail  contour. 

The  nozzle-base-boattail  assembly  is  first  modeled  as  a two  inde- 
pendent variable  problem  by  fixing  the  width  of  the  base  region  to  be 
the  minimum  allowed  width  acceptable  in  the  final  design  (which  may  be 
zero  if  a sharp  trailing  edge  is  permissible).  Consequently,  specifi- 
cation of  the  nozzle  exit  lip  radius  y^-.^  fixes  the  value  of  the  boat- 
tail  exit  radius  yeb^  so  that  the  boattail  contour  is  completely  spec- 
ified. The  nozzle  throat  attachment  angle  0^^  and  exit  lip  radius  yen 
remain  free  to  vary.  This  constitutes  an  optimization  in  two  indepen- 
dent parameters,  which  is  carried  out  until  convergence  is  attained. 

This  two-dimensional  optimization  yields  a good  initial  estimate  for 
the  final  nozzl e-base-boattai 1 maximum  thrust  contour. 

The  problem  is  then  redefined  to  be  an  optimization  in  three  inde- 
pendent variables,  with  the  nozzle  and  boattail  exit  radii  no  longer 
related.  The  three-diinensional  optimization  is  then  carried  out  until 
the  maximum  thrust  producing  contour  is  obtained. 

Two  special  cases  may  arise  that  reduce  the  optimization  problem 
from  a three-dimensional  problem  to  either  a two-dimensional  problem 
or  a one-dimensional  problem.  The  first  case  arises  when  an  optimization 
step  requires  the  boattail  exit  radius  ygb  To  become  less  than  the  nozzle 
exit  radius  yen  (or  less  than  yen  + Ay^,  where  Ayp  is  the  minimum  allow- 
able base  width).  In  that  case,  the  boattail  exit  radius  yeb  is  spec- 
ified as  yeb  = 7en  Ayb*  and  the  solution  continues  as  a two-dimensional 
optimization  problem  with  Gan  and  yen  as  the  independent  variables.  The 
second  case  arises  when  an  optimization  step  requires  a nozzle  exit 
radius  Vgp  greater  than  the  maximum  radius  of  the  boattail.  In  that 
case,  yen  (or  yen  + Ayb)  and  yeb  are  set  equal  to  the  maximum  boattail 
radius,  and  the  solution  continues  as  a one-dimensional  optimization 
problem  with  Can  as  the  single  independent  variable.  The  program  can 
efficiently  perform  an  axial  iteration  optimization  in  one-dimension, 
as  encountered  in  a nozzle  optimization  with  fixed  exit  radius,  boattail 
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optiiiii zdtion  with  fixed  exit  radius,  or  nozzle-base-boattail  assembly 
optimization  with  fixed  boattail  contour  and  nozzle  exit  radius. 

The  analysis  program  is  capable  of  calculating  the  thrust  for 
nozzle-base-boattail  assemblies  that  have  crossing  exit  radii.  That 
situation  can  not  occur  physically,  but  may  occur  during  the  calculation 
of  the  function  derivatives  required  to  compute  an  optimization  step. 

A parametric  study  of  the  thrust  as  a function  of  the  wall  geometry 
for  a nozzle-base-boattail  assembly  was  conducted  to  gain  insight  into 
the  design  of  maximum  thrust  nozzle-base-boattail  contours. 

The  study  was  conducted  for  a nozzle  having  an  upstream  radius  of 
curvature  = 3.0  in.,  a downstream  radius  of  curvature  = 0.5  in., 
a throat  radius  ytn  = 1-0  in.,  and  a length  Xen  = 10.0  in.  The  prop- 
erties of  the  gas  flowing  in  the  nozzle  were  specific  heat  ratio  y = 1.2, 
gas  constant  R = 60.0  (ft-lbf)/(lbm-R) , stagnation  temperature  T = 6000  R, 
stagnation  pressure  P = 1000  psia,  and  ambient  pressure  Pamb  ~ O-O-  The 
boattail  had  a length  xeb  = 10.0  in.,  an  initial  radius  y^b  = 9.0  in., 
and  an  initial  expansion  contour  radius  ptd  = 5.0  in.  The  properties 
of  the  air  flowing  around  the  boattail  were  specific  heat  ratio  y = 1.4, 
gas  constant  R = 53.3  (ft-lbf)/(lbin-R) , stagnation  temperature  T = 600  R, 
and  stagnation  pressure  P = 12.0  psia.  Lines  of  constant  thrust  were 
determined  for  the  nozzles  and  the  boattails  (see  Figs.  3 and  4).  It  is 
seen  that  the  thrust  varies  smoothly,  and  that  there  is  one  clearly 
defined  global  maximum  thrust  for  the  nozzle.  A difference  table  re- 
duction of  the  data  showed  that  locally  both  the  nozzle  thrust  and  the 
boattail  thrust  varied  approximately  quadratical ly. 

Nonlinear  programming  optimization  techniques  may  be  divided  into 
two  broad  categories;  dervi ati ve-free  methods  and  derivative  methods. 
Derivative-free  methods  employ  systematic  searches  through  the  domain 
of  the  objective  function,  in  this  case  the  total  thrust  of  the  nozzle- 
base-boattail  assembly.  These  methods  seek  the  extremum  by  comparing 
several  function  values  (i.e.,  the  thrust)  for  various  combinations  of 
the  independent  variables  Yen’  ^eb)’  moving  to  the  point  having 

the  greater  thrust.  The  idea  is  to  step  to  the  absolute  extremum  in  an 
organized  fashion,  without  recalculating  the  objective  function  at  the 
same  point  more  than  once,  and  without  sampling  the  objective  function 
prohibitively  often.  This  class  of  methods  is  most  useful  when  the 
objective  function  is  a function  of  many  variables. 

Derivative  methods  employ  more  infortietion  about  the  function  than 
a simple  comparison  of  the  function  value  at  two  points.  They  use 
derivative  information  about  the  objective  function,  and  are  most  effec- 
tively employed  where  the  objective  function  is  smooth  and  depends  on 
a few  variables.  This  class  of  methods  was  the  more  successful  in  the 
present  investigation. 

The  general  procedure  of  any  optimization  method  is  to  begin  with 
a starting  point  (base  point).  The  algorithm  then  selects  a direction 
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and  a distance  for  moving  the  base  point  toward  the  extremum.  The 
number  of  base  point  moves  and  the  number  of  function  evaluations 
required  for  each  base  point  move  determines  the  efficiency  of  the 
algorithm,  since  a function  evaluation  requires  much  more  computational 
effort  than  the  optimization  logic. 

Inherent  in  all  optimization  procedures  is  a weakness  at  abrupt 
folds  of  the  objective  function  surface  and  along  ridges  and  plateaus. 
Fortunately,  the  objective  function  considered  in  the  present  investi- 
gation (i.e.,  the  thrust)  does  not  have  any  of  these  features,  thus 
simplifying  the  algorithms. 

Each  optimization  method  requires  an  initial  estimate  of  the 
maximum  thrust  contour.  The  method  then  systematically  revises  the 
estimate,  moving  to  the  extremum  of  the  objective  function.  This 
initial  estimate  should  be  reasonably  close  to  the  final  solution,  so 
that  the  optimization  does  not  require  a prohibitive  number  of  steps, 
or  the  analysis  does  not  encounter  flow  separation  or  strong  internal 
oblique  shock  waves.  One-dimensional  gas  dynamic  relations  for  the 
area  ratio  required  to  expand  the  nozzle  flow  to  the  freestream  pressure 
are  sufficient  to  estimate  the  nozzle  exit  radius.  Some  suitably  small 
annular  base  width  specifies  the  initial  boattail  exit  radius. 

'■  Finite  perturbations  for  each  of  the  independent  parameters  are 
used  in  forming  the  finite-difference  derivative  approximations  used 
in  each  of  the  optimization  methods.  These  perturbations  must  be  large 
enough  so  that  the  corresponding  thrust  increments  are  not  influenced 
by  the  accuracy  of  the  numerical  algorithm,  yet  small  enough  to  accu- 
rately approximate  local  derivatives  of  the  thrust.  Also,  each  per- 
turbation should  change  the  thrust  by  roughly  the  same  amount.  Typical 
values  for  the  perturbations  are  2.0  deg  for  the  nozzle  attachment 
angle,  and  20  percent  of  the  nozzle  throat  radius  for  the  nozzle  and 
boattail  exit  radii. 


A . Axial  Iteration 

Axial  iteration  is  a multidimensional  line  search  algorithm.  It 
treats  the  objective  function  as  a function  of  several  coordinates,  then 
performs  a line  search  in  each  of  the  coordinate  directions  successively. 
A line  search  finds  the  distance  for  a base  point  move,  given  a search 
direction.  This  may  be  envisioned  as  passing  a plane  through  the 
objective  function  surface,  through  the  base  point,  and  along  the  search 
direction.  The  line  search  seeks  to  move  to  the  maximum  on  this  slice. 
Given  a priori  that  the  surface  is  approximately  quadratic  locally,  the 
line  search  perturbs  the  independent  variable  along  the  search  direction 
by  an  equal  increment  on  both  sides  of  the  base  point.  First  and  second- 
derivative  information  is  formed  using  these  function  values  for  base 
point  moves  in  one-dimension  (see  Newton's  method),  the  dimension  suc- 
cessively being  each  of  the  independent  parameters.  The  new  base  point 
estimate  is  then  found,  and  the  objective  function  is  evaluated  at  that 
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point.  The  base  point  moves  to  tiie  niaximuni  of  the  points  thus  far 
evaluated,  and  proceeds  to  the  next  independent  variable.  Convergence 
is  attained  when  a base  point  niove  ciianges  the  objective  function  less 
than  a specified  relative  tolerance. 

Axial  iteration  requires  3n  (where  n is  the  number  of  independent 
parameters)  function  evaluations  per  pass  through  the  algorithm,  and 
does  not  have  difficulty  with  poor  initial  estimates  for  the  maximum 
thrust  contour.  Further,  this  method  has  one  of  the  good  features  of 
the  nonderivative  methods,  in  that  it  steps  to  the  maximum  of  the  points 
evaluated,  not  necessarily  the  predicted  point.  It  is  a very  simple 
method,  and  converges  very  rapidly. 


B . Method  of  Steepest  Descent 

Method  of  steepest  descent,  or  hill  climbing,  is  a first  derivative 
method.  By  perturbing  the  objective  function  n times  (where  n is  the 
number  of  independent  variables)  about  the  base  point  (n  + 1 function 
evaluations),  a finite  difference  approximation  to  the  gradient  may  be 
calculated.  The  search  direction  is  then  specified  along  the  gradient 
direction.  The  algorithm  liiarches  stepwise  in  that  direction,  seeking 
the  extremum,  and  moves  the  base  point  to  that  extremum. 

If  the  first  step  after  computing  an  approximation  to  the  gradient 
does  not  increase  the  objective  function,  which  normally  occurs  near 
the  maximum,  the  perturbation  sizes  and  the  step  size  are  halved.  A 
new  approximation  to  the  gradient  is  then  computed.  The  procedure  is 
repeated  from  successive  base  points  until  convergence  is  achieved. 

This  method  is  relatively  simple,  and  convergence  is  monotonic. 

The  principal  drawback  to  the  method  of  steepest  descent  is  that 
the  search  can  sometimes  zig-zag  back  and  forth  in  a so-called  "hem- 
stitching pattern"  (9).  This  is  especially  true  along  ridges  in  the 
objective  function  surface.  Consequently,  the  method  may  be  very 
inefficient  and  require  a large  number  of  steps.  However,  the  algorithm 
will  eventually  converge  to  the  extremum. 

This  method  is  reconriended  primarily  in  cases  where  a good  initial 
estimate  of  the  wall  geometry  is  not  available  or  to  check  the  result 
obtained  by  any  other  method  to  ensure  that  the  absolute  maximum  is 
attained.  After  several  base  point  moves,  one  of  the  other,  faster 
converging,  methods  may  be  employed  to  complete  the  solution. 


C.  Newton's  Method 

Newton's  method  is  a second-derivative  method.  The  objective 
function  is  approximated  by  a Taylor  secies  expansion  about  the  base 

rgl  - 
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point  which  is  tr'uncated  after  the  second-derivative  terms.  In  one- 
dimension,  the  following  equation  is  obtained. 


f(x)  = f(a)  + f^(a)(x-a)  + ~r  <"^^(a)(x-a)^ 


(8) 


wtiere  f is  the  objective  function,  a is  the  base  point,  and  the  deriva- 
tives are  evaluated  at  the  base  point,  a.  Taking  the  derivative  of 
equation  (8)  with  respect  to  x and  noting  that  f(a),  f (a),  and  f (a) 
are  constant  yields  ^ 

At  a stationary  point  x = x*  (maximum,  minimum,  or  saddle  point), 
is  identically  zero.  Consequently,  at  a stationary  point,  equation^(9) 
becomes 

0 = fx(a)  + no) 

Solving  equation  (10)  for  x*  yields 

X*  = a - f^(a)/f^^(a)  (11) 

If  f(x)  is  not  quadratic,  x*  is  not  exactly  the  extremum.  Usually, 
however,  it  is  an  improvement,  and  convergence  is  approached  iteratively. 

In  the  case  of  more  than  one  dimension,  x*,  a,  and  fx(a)  are 
replaced  by  column  vectors,  and  fxx(a)  is  an  n x n matrix  of  second 
derivatives,  termed  the  Hessian.  Thus,  equation  (11)  becomes 

X*  = a - (12) 

Solving  for  x*  in  equation  (12)  requires  inversion  of  the  Hessian, 
which,  in  general,  may  present  a numerical  accuracy  problem.  In  the 
present  problem,  where  n ^ 3,  the  Hessian  has  proven  to  be  nonsingular 
and  well  conditioned,  and  no  problem  arises. 

Newton's  method  converges  very  rapidly.  In  fact,  it  converges  in 
exactly  one  pass  for  a quadratic  function.  However,  it  does  require  a 
relatively  large  number  of  function  evaluations  for  each  base  point 
move  [(n2  + 3n  + 2)/2  evaluations/move]. 

When  applying  Newton's  method,  it  may  be  possible  to  employ  the 
Hessian  computed  for  a given  base  point  for  several  base  point  moves 
before  recomputing  second  derivatives.  In  the  case  of  a nearly  quadratic 
function,  the  second  derivatives  are  approximately  constant,  and  the 
procedure  may  be  useful.  The  thrust  function  considered  in  the  present 
study  was  not  close  enough  to  a quadratic  function  globally,  and  the 
method  did  not  converge  unless  second  derivatives  were  computed  at  each 
step. 


14 


i 


Furthet',  when  the  objective  function  is  not  ylobully  quadratic, 
Newton's  method  may  predict  bad  base  point  moves,  particularly  if  the 
initial  estimate  is  far  from  the  extremum.  Hence,  there  is  special 
incentive  for  a good  firsc  guess. 

Newton's  method  perforins  quite  well  in  the  present  problem.  Bard 
(9)  states  that  it  is  relatively  common  to  find  cases  where  Newton's 
method  may  be  Z5,000  times  more  efficient  than  the  method  of  steepest 
descenti  Though  not  outperforming  the  other  methods  nearly  so  dramati- 
cally in  the  present  study,  it  is  clearly  an  efficient  method. 


15 


SECTION  III 


RESULTS 


The  optimization  methods  presented  in  Section  II  have  been  applied 
to  the  design  of  maximum  thrust  nozzle  contours,  boattail  contours,  and 
nozzle-base-boattail  contours.  These  contours  are  modeled  as  first- 
and  second-order  polynomials.  The  nozzle  thrust  predicted  by  the  maximum 
thrust-producing  second-order  polynomial  is  comparable  to  that  predicted 
by  the  calculus  of  variations,  which  places  no  restrictions  on  the 
functional  form  of  the  wall  contour.  A study  comparing  the  results  of 
the  two  approaches  was  conducted  for  different  length  nozzles  operating 
at  two  ambient  pressures. 

A method  of  nozzle  optimization  employing  the  calculus  of  variations 
is  based  on  forming  a functional  consisting  of  the  following  terms: 

(1)  the  term  to  be  maximized  (the  axial  pressure  forces  along 
the  contour  y = y(x) ) , 

(2)  the  design  constraint  along  the  nozzle  wall  y = (such 

as  constant  nozzle  length  or  surface  area), 

(3)  the  constraint  that  the  contour  y = y (x)  is  a streamline, 

and  ^ 

(4)  the  constraint  that  the  governing  partial  differential  equations 
must  be  satisfied  in  the  region  of  the  flow  field  influenced  by 
the  nozzle  contour. 

Each  of  the  constraints  is  multiplied  by  a Lagrange  multiplier,  and  a 

sum  of  terms  (1)  through  (4)  is  formed.  No  constraint  is  placed  on 

the  form  of  y = y,  ,(x)  • 
w 

Following  the  formalism  of  the  calculus  of  variations,  the  first 
variation  of  the  functional  is  determined  and  set  equal  to  zero.  A set 
of  conditions,  known  as  Euler  equations,  transversal ity  equations,  and 
corner  conditions,  results  from  setting  the  first  variation  equal  to 
zero.  These  design  equations  guarantee  a maximum  thrust  contour  for 
the  assumed  constraints.  Because  the  design  equations  are  nonlinear,  an 
iterative  solution  is  performed.  A trial  contour,  y = , is  assumed, 

the  flow  field  is  determined  by  the  method  of  characteristics,  and  the 
design  equations  are  solved  for  the  Lagrange  multipliers.  During  the 
above  procedure,  one  transversal i ty  equation  along  the  nozzle  wall  is 
not  employed.  That  design  equation,  expressed  in  terms  of  the  flow 
field  properties,  the  nozzle  contour  y = yy^(x)»  ^nd  the  Lagrange  multi- 
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pliers,  is  then  checked  to  see  it  it  is  satisfied.  If  it  is  satisfied, 
the  assumed  contour  is  the  maximum  thrust  no/zle  contour.  If  it  is  not 
satisfied,  the  contour  is  iiKidiried  and  the  entire  procedure  is  repeated 
for  the  modified  contour.  Inis  relaxation  of  the  nozzle  contour  is 
repeated  until  the  design  condition  is  satisfied  within  a specified 
tolerance.  Typically  five  to  ten  modifications  of  the  nozzle  contour  are 
required.  For  each  modification,  the  nozzle  flow  field  and  Lagrange 
multiplier  field  must  be  calculated;  the  Lagrange  multiplier  field  calcu- 
lation requiring  approximately  60  percent  as  much  computation  time  as 
that  required  for  the  nozzle  flow  field  calculation.  Thus,  the  indirect 
approach  requires  from  eight  to  sixteen  times  the  computation  time 
required  for  one  nozzle  flow  field  calculation. 

The  calculus  of  variations  is  an  indirect  approach  to  finding  maxi- 
mum thrust  nozzle,  boattail,  and  nozzle-base-boattai 1 contours.  The 
advantage  of  the  indirect  approach  is  that  no  restrictions  are  placed  on 
the  general  forms  of  the  contours  y = yw(x).  The  contours  are  free  to 
assume  any  shape  required  to  yield  maximum  thrust.  The  disadvantages  are 
that  a tremendous  amount  of  time  is  required  to  formulate  and  solve  the 
variational  problem  to  obtain  the  design  equations,  and  to  develop  a 
computer  program  to  solve  for  the  Lagrange  multiplier  field  for  a maxi- 
mum thrust  nozzle.  This  problem  is  further  complicated  by  consideration 
of  base  and  boattail  effects.  Even  if  these  could  be  accounted  for,  a 
subject  of  some  speculation,  the  entire  problem  would  need  to  be  entirely 
reworked  if  changes  were  made  to  the  gas  dynamic  model,  such  as  the 
inclusion  of  boundary  layer  drag  or  nonequilibrium  chemistry  effects. 

In  spite  of  these  disadvantages,  the  indirect  method  is  of  great  value 
in  that  it  yields  absolutely  the  best  solution  to  the  specified  problem. 

The  concept  of  modeling  the  contours  of  a complete  nozzle-base- 
boattail  assembly  as  polynomials,  and  performing  a direct  optimization 
to  find  maximum  thrust  has  proven  to  be  practical,  and  an  example  is 
presented  to  demonstrate  the  application  of  the  optimization  techniques. 


1 . Nozzle  Design 

Sample  nozzle  optimizations  were  carried  out  using  the  three 
methods  outlined  in  Section  II  for  the  nozzle  illustrated  in  Fig.  3. 
Their  paths  to  the  extremum  is  shown  in  Fig.  5,  which  illustrates  that 
each  of  the  techniques  has  the  following  properties. 

a.  Each  converges. 

b.  Each  converges  to  the  same  stationary  point. 

c.  The  stationary  point  converged  upon  is  the  global  maximum. 

d.  Each  method  converges  to  the  global  maximum  efficiently, 
without  an  excessive  number  of  base  point  moves. 
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NOZZLE  ATTACHMENT  ANGLE  DEGREES 


Thus,  several  practical  optiinization  tcctm  icaes  d»'e  available  to  find 
the  nozzle  wall  cieonietry  tfiai  vie''ds  iiiaxii.ia;a  hru..t. 

To  demonstrate  that  -I'Jddratic  wall  contours  develop  thrust 
comparable  to  that  developed  by  the  calculus  of  variations  contours, 
a parametric  study  was  performed  for  a variety  of  nozzle  lengths  and 
two  ambient  pressures.  The  results  of  this  study  are  presented  in 
Table  1.  The  calculus  of  variations  procedure  einployed  in  the  comparison 
is  that  developed  by  Rao  (10). 

The  nozzle  parameters  for  the  study  are  upstream  radius  of  curvature 
wtij  = 3.0  in.,  downstream  radius  of  curvature  = 0.5  in.,  throat  radius 
yi;  = 1.0  in.,  specific  heat  ratio  i - 1.3,  gas  constant  k ^ 60  (ft-lbf)/ 
(Ibm-R),  ftacituUion  temperature  “T  - 6000  R,  and  stagnation  pressure 
P = 1000  psia. 

The  tin  usts  delivered  by  the  contours  designed  by  tne  calculus  of 
variationsare  comiiared  wits  the  <hrusts  delivered  by  the  contours  ae- 
signed  by  the  direct  optimization  of  the  second-order  polynomial  nozzle 
wall.  Both  methods  predicted  essentially  the  same  niaxinium  obtainable 
thrust,  with  the  direct  optimization  -'lethod  not  having  the  analytical 
difficulties  of  vr,e  calculus  of  variations  approacn.  Slight  variations 
of  the  maximum  ttirust  predicted  by  trie  two  approaches  are  acceptable 
because  they  use  different,  aithough  similar,  flow  field  analysis 
programs,  and  the  direct  optimizations  were  terminated  by  nonzero  relative 
tolerances . 

Overall  ayreci'ent  becv.een  the  two  approaches,  that  is,  the  calculus 
of  variations  anc  the  direct  optiisization  of  a second- order  polynomial 
wall  nozzle,  is  very  good.  For  the  case  of  zero  ambient  pressure,  agree- 
ment for  all  cases  is  within  0.2  percent,  which  is  the  approximate 
numerical  resolution  of  the  method  of  cnaracteristic  flow  field  calcula- 
tion algorithm,  fne  results  for  an  ambient  pressure  of  5.0  psia  are  less 
compatible,  in  most  part  due  to  the  second-order  polynomial  nozzle  wall 
approximation  becoming  less  good.  Thus  second-order  polynomials  effec- 
tively approximate  up  .i mum  thrust  nozzle  contours,  ju'^tifying  this 
approach. 


2.  iiq_zz2e-Base-^^oaitdi  1 

The  program  is  capable  of  conou^’ing  the  design  of  a nozzle-base- 
oodttdi!  assei'bly  whicn  yields  maximum  tnrusi.  This  is  of  great  practical 
value  in  that  it  dictates  the  geometry  which  can  produce  maximum  per- 
turmance,  without  necessi cati iig  excessive  experience  or  experimental 
testing.  The  optimization  nodels  the  assembly  as  having  a second-order 
polynomial  nozzle  contour,  a first-order  polynomial  (conical)  boattail 
contour,  and  an  annular  tiase.  This  geometry  is  uniquely  specified 
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TABLL  1.  COMPARISON  OF  THE  PERFORMANCE  OF  :iAXIMUM  THRUST 
NOZZLE  CONTOURS 


Thrust,  Ibf 


Nozzle 

Length 

(in.) 

Amb i en  t 
Pressure 
(psid) 

Cal cu 1 us 
of  Variations 
Prediction 

Second-Order 
Polynomial  Wall 

AF,  percent 

1 .5475 

G.O 

4556.0 

4555.7 

0.01 

2.3215 

0.0 

4802 . 7 

4798.1 

0.10 

3.0706 

0.0 

4981.4 

4979.0 

0.05 

3.8205 

0.0 

5124.5 

5120.9 

0.07 

4.5731 

0.0 

5241  .7 

5234.9 

0.13 

6.2965 

0.0 

5433.1 

5428.7 

0.08 

7.9759 

0.0 

5569.6 

5561 .2 

0.15 

9.6527 

0.0 

5669.4 

5662.5 

0.12 

12.746 

0.0 

5806.7 

5801 .2 

0.09 

16.9164 

0.0 

5935.7 

5926.6 

0.15 

20.5326 

0.0 

6014.4 

6003.2 

0.19 

24.4965 

0.0 

6081.4 

6068.2 

0.21 

1.4139 

5.0 

4468.4 

4453.1 

0.34 

2.8861 

5.0 

4880.4 

4848.3 

0.66 

7.9850 

5.0 

5311.8 

5300.1 

0.22 

11.836 

5.0 

5379.0 

5367.4 

0.22 

16.060 

5.0 

5394.9 

5377.6 

0.32 
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knowing  nozzle  and  boattail  lengths,  initi'l  expansion  contours,  nozzle 
attdchiiient  angle  nozzle  exit  radius  ypn-  and  boattail  exit  radius 

yeb-  The  last  three  parameter:,  are  varied  in  sucii  a way  as  to  find  that 
unique  combination  ( an.  aen’  Veb)  which  yields  maximum  thrust. 

There  are  several  eonstraints  on  these  variables,  besides  the 
minimum  and  maximum  values  for  each  necessary  to  bracket  the  extremum. 

The  nozzle  and  boattail  exits  must  not  cross,  a physical  impossibility. 
This  has  been  taken  one  step  further  in  that  they  can  be  specified  to 
be  at  least  some  specified  minimum  base  width  (Ay5)  away  from  each  other. 
If  the  optimization  violates  this,  the  two  exit  radii  become  a fixed 
distance  apart,  and  the  optimization  proceeds  with  one  less  independent 
variable  (n  ^ 2).  The  nozzle  exit  radius  may  also  be  constrained  to 
be  loss  than  some  maximum  value.  Tfiis  could  be  encountered  if  the 
nozzle  exit  radius  became  larger  than  that  of  the  airframe.  This  causes 
tne  nozzle  exit  radius  to  become  fixed  and  the  optimization  proceeds  with 
one  less  independent  variable  (n  = 1). 

A nozzle-base-boatta i 1 optimization  is  presented  in  Table  2 to 
illustrate  the  application  of  the  program.  Each  of  the  three  methods 
in  Section  II  is  applied  to  the  same  case.  The  nozzle  is  the  same  as 
illustrated  in  Fig.  3.  The  boattail  is  the  same  as  illustrated  in  Fig. 

4 except  that  it  tias  an  outer  radius  - 5.3  in.  Note  that  axial 
iteration  and  Newton's  method  converge  comparably,  while  the  method  of 
steepest  descent  converges  to  a larger  tolerance  much  more  slowly. 
Newton's  method  is  more  prone  to  failure  than  axial  iteration,  because 
it  requires  a better  estimate  of  the  partial  derivatives  to  take  a 
successful  step. 


Thus,  maximum  thrust  nozzle-hase-boattail  contours  may  be  approxi- 
mated as  polynomials  and  may  be  directly  optimized  efficiently  using 
several  techniques.  The  order  of  the  approximating  polynomials  seems 
to  be  effective.  The  nozzle  contour  has  been  shown  conclusively  to  be 
well  approximated  by  a second-order  polynomial.  Since  the  boattail 
typically  contributes  an  order  of  magnitude  less  thrust  than  does  the 
nozzle,  it  was  modeled  as  conical.  Thus,  there  are  as  many  as  three 
independent  parameters,  since  the  attachment  angle  may  be  inferred  from 


7ABLE  2.  RESULTS  OF  A MOZZLE-BASE-BOATTAIL  Tlll’.EE-jlMEriSIONAL  OPTIMIZATION 


Function 

0 , 

an 

^en’ 

9 u. 
ab 

■^eb’ 

F,  Ibf 

Method 

Step 

n 

Evaluations 

deg 

in. 

deg 

in. 

Axial 

0 

2 

1 

30.000 

4.00 

7.059 

4.100 

5710.71 

iteration 

1 

2 

7 

27.223 

4.296 

5.344 

4.396 

5734.72 

2 

2 

13 

29.365 

4.377 

4.906 

4.477 

5737.39 

3 

2 

19 

29.579 

4.391 

4.722 

4.491 

5737.46 

0 

3 

20 

29.578 

4.391 

4.130 

4.591 

5737.65 

1 

3 

29 

28.700 

4.369 

3.780 

4.650 

5737.74 

Steepest 

0 

2 

1 

30.000 

4.000 

7.059 

4.100 

5710.71 

descent 

1 

2 

9 

29.947 

4.497 

4.094 

4.597 

5735.93 

2 

2 

12 

29.947 

4.497 

4.094 

4.597 

5735.93 

3 

2 

16 

29.938 

4.448 

4.385 

4.548 

5736.09 

4 

2 

19 

29.932 

4.448 

4.395 

4.548 

5736.09 

5 

2 

24 

29.888 

4.443 

4.413 

4.543 

5736.70 

0 

3 

25 

29.888 

4.443 

3.923 

4.643 

5736.27 

1 

3 

30 

29.888 

4.443 

3.823 

4.643 

5736.27 

2 

3 

34 

29.888 

4.443 

3.823 

4.643 

5736.27 

3 

3 

39 

29.973 

4.427 

3.889 

4.632 

5736.29 

Newton 's 
method 

0 

2 

1 

30.000 

4.0 

7.059 

4.100 

5710.71 

1 

2 

7 

29.597 

4.235 

5.653 

4.335 

5734.89 

2 

2 

13 

29.552 

4.349 

4.974 

4.449 

5737.30 

3 

2 

19 

28.710 

4.399 

4.674 

4.499 

5737.47 

0 

3 

20 

28.710 

4.399 

4.082 

4.599 

5737.64 

1 

3 

30 

28.551 

4.347 

3.661 

4.670 

5737.76 

SUCTION  IV 


CO'."  ' "'JONS 


A method  has  been  develoi^ed  tor  the  design  of  inaximuin  thrust  nozzles, 
boattdils,  and  nozzle-base-boattai 1 contours.  Triese  contours  may  be 
quickly  and  efficiently  found  using  any  of  several  direct  optimization 
techniques,  providing  a practical  anu  dependable  basis  for  design. 

Nozzle  and  boattail  op ti fi;i z.n  ion:-  seieci  second-order  polynomial 
wall  contours.  Nozz  Ic-base-boatta  i opMmizations  assuiiie  a first-order 
polynomial  (conical)  boatta i I contojt  and  a second-order  polynomial 
nozzle  contour.  Better  approximations  for  botli  nozzle  and  boattail 
maximum  thrust  contours  may  be  oiitained  using  higher-order  polynomial 
walls,  and  hence  degrees  of  freedom,  fach  additional  degree  of  freedom 
allows  a better  fit  of  the  general  surface  which  yields  maximum  thrust. 
However,  each  additional  degree  of  freedom  also  increases  the  number  of 
independent  contour  parameters  required  to  uniquely  specify  the  contour 
geometry,  geometrically  increasing  the  number  of  function  evaluations 
necessary  to  perform  an  optimization.  The  second-order  polynomial  wall 
nozzle  and  conical  wall  boattail  were  selected  as  an  effective  compromise. 
The  polynomials  are  shown  to  yield  substantial ly  identical  thrust  to  that 
produced  by  maximum  thrust  contours  designed  by  the  calculus  of  varia- 
tions for  unconstrai ned  wal ! geometry. 

Three  techniques  were  oevelopeii  for  determining  the  wall  polynomial 
coefficients.  The  first  method  is  an  axial  iteration.  It  is  a fast 
converging  method,  and  has  tne  least  difficulty  with  a poor  initial 
estimate  for  the  solutions  of  the  methods  studied.  This  piethod  also 
can  take  steps  from  simple  sampling,  if  terms  of  higher  than  second 
order  dominate  the  objective  function  and  second  and  first  derivative 
information  yields  a poor  step. 

The  second  technique  is  the  method  of  steepest  descent.  This 
method  converges  monotonical ly  to  the  extremum  and  cannot  overstep  it. 

It  may,  however,  require  an  exhurb inant  number  of  steps  to  attain  the 
maximum.  It  can  step  directly  away  from  a poor  initial  estimate  of  the 
maximum,  or  converge  to  the  extreM'.um  frcm  an  approximate  solution  ob- 
tained from  one  of  the  otiier  methods.  This  algorithm  is  primarily 
intended  as  a classical  reference. 

The  third  method  is  Newton's  method.  This,  too,  is  one  of  the 
classical  optimization  methods,  and  woi'ked  very  well  in  the  present 
study.  Newton's  method  does  reouire  accurate  derivative  approximations, 
and  care  must  be  taken  to  be  sure  finite-difference  increments  are  large 
enough  to  be  beyond  the  numerinal  resolution  of  the  analysis,  yet  small 
enough  to  detect  local  functicna''  trends. 


All  three  optimization  techniques  efficiently  select  the  wall 
geometry  that  produces  maximum  thrust  in  a nozzle-base-boattai 1 assembly. 
This  makes  for  a simple  and  practical  approach  to  the  design  of  maximum 
thrust  nozzle-base-boattai 1 contours. 
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APPcfJDICES 


APrCNDIX  A.  CHARAriFRISTIC  AND  COMPATIBILITY  EQUATIONS 
The  following]  equations  qovern  steady  two-dimensional  isentropic 


(i.e.,  rotational)  flow, 

ou^  + , Vy  + v/y  - 0 (A-1) 

puu^  t-  vUy  = 0 (^"2) 

ruv^  + pvVy  + = 0 (A-3) 

aPx  vpy  - a^upx  - a^vp^  = 0 (A-4) 


where  6=0  for  planar  flows  and  6 = 1 for  axisynmetric  flows.  For 
isentropic  flow,  the  speed  of  sound,  a,  is  a known  function  of  the 
static  pressure  p and  the  density  p.  Thus, 

a = a(p,p)  (A-5) 

Equations  A-1  througii  A-4  comprise  a system  of  four  quasi-linear 
nonhomogeneous  first-order  partial  differential  equations  for  the  four 
variables  u,  v,  p,  and  p. 

This  system  of  partial  differential  equations  may  be  replaced  by 
an  equivalent  system  of  four  compatibility  equations  which  are  valid 
along  three  distinct  characteristic  curves:  the  streamline  and  the  two 

Mach  lines.  The  compatibility  equations  are  total  differential  equa- 
tions, as  are  the  equations  defining  the  characteristic  curves.  Thus, 
a much  simpler  numerical  method  can  be  employed  to  solve  the  system  of 
total  differential  equations. 

The  characteri s tic  and  compatibil i ty  equations  for  the  above 
system  of  equations  are  found  by  multiplying  equations  A-1  through  A-4 
by  the  unknown  parameters  o]  through  .'4,  respectively,  and  summing. 

Thus , 


^(A-1)  + -,^(A-2)  + -,3(A-3)  + >^(A-4)  = 


.’•J 


0 


(A-6) 


After  this  sum  is  formed,  the  coefficients  of  the  x derivatives  of  u, 
V,  p,  and  p are  factored  out  to  yield 


pvo. 


(pOi  + pUOo)  U + ------ 

1 ^ 2M  X pa,  + p 


po^  + PUO2  y 


+ (puo^) 


''x" 


pa-|  + pva^ 
pua. 


(a^  + uo^) 


03  + VO, 

P + P 

^x  02  + ua^  '^y 


(uo^  - a UO4) 


P + 

'^x 


va-|  - a va^ 

—~2  ^^y 

uo-j  - a uo^ 


+ a^6ov/y  = 0 (A-7) 


The  slopes  of  the  characteristic  curves,  dy/dx  = A,  are  then  the 
coefficients  of  the  y derivatives  of  u,  v,  p,  and  p in  equation  A-7. 
Thus , 


VO, 


A = 


a,|  + uo^ 


a-|  + V03  03  + vo^  vo,|  - a vo^ 

UO-,  Oo  + uo.  2 

3 2 4 uo-|  - a uo^ 


(A-8) 


If  u,  V,  p,  and  p are  assumed  to  be  continuous  functions,  then 
du/dx  = Ux  + 'Uy,  etc.,  and  the  terms  in  brackets  in  equation  A-7  may 
be  replaced  by  their  equivalent  values,  du/dx,  etc.  Thus,  equation 
A-7  becomes 


p(o,|  + uo2)du  + puo3dv  + (02  + uo^)dp  + 

u(a,|  - a^o^)dp  + o,|(6pv/y)dx  = 0 (A-9) 

Equation  A-9  is  the  compatibility  equation,  which  is  valid  along  the 
characteristics  defined  by  Equation  A-8.  It  remains  to  eliminate  a, 
to  o^  from  equation  A-8  and  A-9. 

Equation  A-8,  when  written  as  a system  of  equations  considering 
a-|  through  0^  as  the  unknowns,  becomes 


Oi(A)  + 02(uA  - v)  + 03(0)  + 0^(0)  = 0 

(A-10) 

Pl(-1)  + 02^°^  03(uA  - v)  + 0^(0)  = 0 

(A-11) 

o,j(0)  + 02(A)  + 03(-l)  + 04(uA  - v)  = 0 

(A-12) 
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(A-13) 


0^ ( UA  - v)  + 


a^(0)  + 0^(0)  t o^[-a^{uA  - v)]  = 0 


For  equations  A-10  throu'jh  A-13  to  have  any  solution  for  the  o's  other 
than  the  trivial  solution  of  zero,  the  determinant  of  the  coefficient 
matrix  of  equations  A-10  to  A-13  must  be  zero.  Defining  S = (u>  - v), 
that  condition  yields 


-1 

0 

S 


S 0 0 

0 S 0 

A -1  S 

0 0 -a^S 


= 0 


(A-14) 


Expanding  the  determinant  gives 

S^[S^  - a^(l  + A^)]  = 0 


(A-15) 


Equation  A-15  is  a fourth-order  polynomial  equation  in  terms  of  A; 
hence  four  roots,  and  thus  four  characteristic  directions,  should  be 
found.  Two  roots  are  obtained  by  setting  the  leading  term,  S^,  equal 
to  zero.  Thus, 

Ao  = (^)  = ■“  (repeated  twice)  (A- 16) 

0 

which  is  a repeated  root,  and  is  the  differential  equation  for  the 
streamline.  The  subscript  o is  used  to  denote  the  streamline.  Hence, 
the  streamline  is  a dual  characteristic  in  rotational  flow.  The  remain- 
ing two  roots  are  obtained  from  the  quadratic  term  in  equation  A-15, 
which,  when  substituting  S = (uA  - v),  becomes 

(u^  - a^)A^  - 2uvA  + (v^  - a^)  = 0 (A-17) 

Solving  equation  A-17  for  A gives 


A 


t 


uv 


(A-18) 


where  the  subscripts  i correspond  to  the  root  of  the  quadratic  equa- 
tion. Equation  A-18  can  be  simplified  by  making  the  substitution 
u = Vcose,  V = VsinO,  0 = tan"  (v/u)  and  a = sin"l(l/M).  The  result 
is 

A,  = = tan(0  * a)  (A-19) 
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Equation  A-19  is  the  differential  equation  of  the  Mach  lines  for  super- 
sonic flow;  hence,  the  remaining  two  characteristics  for  rotational  flow 
are  the  Mach  lines.  The  + denotes  the  left-running  Mach  line  and  the  - 
denotes  the  right-running  Mach  line. 

Thus,  three  distinct  characteristics  exist  through  each  point  in  a 
rotational  flow:  the  streamline  and  the  two  Mach  lines. 

The  compatibility  equation  valid  on  each  characteristic  is  obtained 
from  Equation  A-9  by  solving  equations  A-10  through  A-13  for  o]  through 
o^.  Equations  A-10  through  A-13  may  be  rewritten  as  follows. 


02  = -(X/S)o^ 

(A-20) 

Oi  = S03 

(A-21) 

(A-22) 

a^So^  = So.| 

(A-23) 

Along  streamlines,  S = (uA  - v)  = 0,  and  equations  A-20  through 
A-23  become 


o 

II 

D 

(A-24) 

O 

II 

(A-25) 

^2^  - 03  = 0 

(A-26) 

0=0 

(A-27) 

Thus,  along  streamlines. 

a-|  = 0 and  = Ao2 

(A-28) 

and  02  and  04  are  unspecified,  and  hence  arbitrary, 
equation  A-28  into  equation  A-9  yields 

Substi tuting 

2 

02[pu  du  + pv  dv  + dp]  + o^[u  dp  - a u dp] 

= 0 (A-29) 

Since  02  and  04  are  both  arbitrary,  their  coefficients  must  be  identi- 
cally zero.  Thus,  along  streamlines 

pu  du  + pv  dv  + dp  = 0 

(A- 30) 

dp  - a^  dp  = 0 

(A-31) 
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On  the  Mach  lines,  the  quadratic  factor  in  equation  A-15  is  zero, 
which  gives 

- a^(1  + X^)  = 0 (A-32) 

2 

Solving  for  (S/a  ) from  equation  A-32  and  substituting  the  result  into 
equation  A-23  yields  equation  A-22,  which  shows  that  only  one  of  these 
equations  is  independent  on  Mach  lines.  Thus,  along  Mach  lines, 

2 

a-|  = Sa^,  and  = (S/a  )a^  (A-33) 

and  o,  remains  arbitrary.  Substituting  equation  A-33  into  equation 
A-9  and  dividing  by  03,  which  is  nonzero  and  arbitrary,  yields  the 
compatibility  equation  valid  on  Mach  lines.  Thus, 

2 

(pv)du^  - (pu)dv^  + [X^  - u(uA^  - v)/a  ]dp^  - 


[v(u>.  - v)/y]dx^  = 0 


(A-34) 


The  subscript  ^ indicates  that  equation  A-34  is  valid  on  left-running 
Mach  lines,  and  the  subscript  - denotes  that  equation  A-34  is  valid  on 
right-running  Mach  lines. 


An  alternate  and  more  useful  form  of  equation  A-34  is  obtained  by 
making  the  substitutions  u = Vcos6,  v = VsinG,  0 = tan“l(v/u),  and 
a = sin-1 (i/M) . Thus  , 


J M^-1 

PV^ 


dp. 


'^yVMcosO  ± a)’ 


dx,  = 0 


(A-35) 


where  the  upper  subscripts  on  dp,  d6,  and  dx  correspond  to  the  upper 
signs  in  the  terms  + d0  and  cos(6  ± a),  and  vice  versa.  Equation  A- 30 
along  a streamline  becomes 

pV  dV  + dp  = 0 (A-36) 

which  is  Bernoulli's  equation. 

Thus,  three  distinct  characteristics  are  obtained,  the  streamline 
and  the  two  Mach  lines,  with  two  cwipati bi 1 i ty  equations  valid  on  the 
streamline  and  one  compatibility  equation  valid  on  each  Mach  line,  a 
total  of  four  compatibility  equations.  This  system  of  characteristic 
and  compatibility  equations  is  sufficient  to  replace  the  original 
system  of  four  partial  differential  equations. 

In  the  numerical  implementation  of  equation  A-35  along  the  Mach 
lines,  left-running  Mach  lines  can  become  essentially  vertical.  In  that 
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case,  the  term  dx  /cos(0  + u)  becomes  indeterminate.  However,  from  equa- 
tion A-19, 

dx , dy , 

cos{0  ± a)  sin(0  ilx)  (A-37) 

Substituting  equation  A-37  into  equation  A-35  yields 


dp^  + de,  + 6 


V 

yVMsin(0  + a) 


dy. 


0 


(A-38) 


In  the  numerical  algorithms  developed  during  this  investigation,  equation 
A-35  is  employed  along  right-running  Mach  lines  [that  is,  dy/dx  = tan 
(6  - a)],  and  equation  A-37  is  employed  along  left-running  Mach  lines 
[that  is,  dy/dx  = tan(o  + a)]. 
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APPENDIX  B.  KLIEGEL'S  TRANSONIC  SOLUTION 


Various  techniques  have  been  developed  to  handle  the  transonic 
flow  region  in  convergent-divergent  noz;^les,  from  which  the  method  de- 
veloped by  Kliegel  (8)  was  selected.  Kliegel  employs  an  expansion  in 
inverse  powers  of  the  normalized  throat  wall  radius  of  curvature 
R(R  = P|j/yt)  • His  method  differs  from  that  developed  by  Hall  (9)  by 
the  use  of  a toroidal  coordinate  system  instead  of  a cylindrical  co- 
ordinate system.  A toroidal  coordinate  system  results  in  an  expansion 
in  1/(R  + 1)  instead  of  1/R  for  the  axial  and  radial  velocity  components. 
This  expansion  allows  for  the  analysis  of  nozzles  having  small  normalized 
throat  wall  radii  of  curvature,  R < 1.  The  solution  derived  by  Kliegel 
is  expressed  in  the  following  form. 

u-i  ( r , z)  1 j-  , 

“ “ ' ♦ (ir”iT " (rT'i7  * «2(t-,z)J+ 


1^ 

(R  + 1)' 


V = 


u^(r,z)  + U2(r,z)  + u^Cr.z)  +•••  (B-1) 


2(R  + 1)  J R + ■>  (R  + 1 


j v^ (r,z)  + V2(r,z) 


1 


(R  4 1)' 


— v^(r,z)  4 I V2(r,z) 


(B-2) 


where 


2(R  4 1)' 
Y 4 1 J 


1/2 


x/yt 


(B-3) 


r = y/y^  (B-4) 

In  Eq.  B-3  z is  the  transformed  normalized  axial  coordinate.  In  Eq.  B-4 
r is  the  normalized  radial  coordinate.  The  complete  third-order  axisy- 
mmetric  solution  derived  by  Hall  and  modified  by  Kliegel  is  as  follows. 
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(B-10) 


V ^ 831  r + L242  1 ^ 2|>2/_t_ 
864  J L 

^ + 75,  + 279  ] 3 [7y  - 3 / 

192  ^ I ^Z 


62.'^  + 51y  + 327  3 


52y^  + 75 V + 279  3 7y  - 3 

192  ^ L 12 

Kliegel's  analysis  is  valid  for  irrotational  flow  of  a perfect 
gas  in  an  axi symmetric  circular  arc  throat. 


APPENDIX  C.  NOMINAL  CASE  AND  CORRELATION 


FACTORS  FOR  BASE  PRESSURE  MODEL 


The  pressure,  p^,,  acting  on  an  annular  base  area  depends  on  the 
boattail  Mach  number  and  static  pressure,  the  nozzle  exit  Mach  number, 
exit  angle  and  static  pressure,  the  specific  heat  ratios  of  the  gases, 
the  stagnation  temperatures  of  the  two  streams,  and  the  ratio  of  nozzle 
exit  radius  to  boattail  exit  radius.  A model  which  accurately  predicts 
the  base  pressure  as  a function  of  all  these  parameters  has  not  been 
developed  yet,  but  the  model  by  Addy  (12)  does  include  the  effects  of 
each  of  these  parameters  and  gives  a reasonable  approximation  to  the 
base  pressure.  Byington  and  Hoffman  (1)  correlated  his  data  by  a curve 
fit  of  the  form 

Ph 

— =f(p  /P  u)n(M  . )n(M  )n(T  ./T  )n(0  ) (C-1) 

Peb  ^en  ^eb  eb  ' en  ' eb  en  en  ' ' 

The  nominal  case  used  for  f(Pen/Peb^  was:  Mg^  = 2.0,  y = 1*2, 

Men  = 2.0,  0en  = 0 deg,  (Teb/Ten)  = 1-0,  and  (yeb/yen)  = 0-6-  This 
resulted  in  the  following  third-order  polynomial  curve  fit. 

T(Pen/Peb)  = 0.27435080  + 0.13825822  x (Pg^/Pg^)  " 

0.01754651  X (Pgp/Pgb)^  + 0.00097156  x (Pg^/Pgb)^  (C-2) 


The  four  n functions  are  correlation  factors  which  are  applied  to 
the  nominal  case.  Each  factor  accounts  for  variations  in  the  parameter 
denoted  by  its  argument  and  was  obtained  by  a curve  fit  of  the  data 
presented  by  Addy  (12).  Correlation  factors  based  on  the  specific  heat 
ratios  and  exit  lip  point  coordinates  have  been  neglected  since  their 
influence  on  the  base  pressure  is  slight.  The  four  n correlation 
factors  are 


n(Mgb)  = 


2.0 


M 


eb 


(C-3) 


3.5  - 2.5/M 

n(M.J  = - 


en' 


M 


(C-4) 


en 


n(0en)  =1-0+ 

The  correlation  factors  n(M  „)  and  n(0„„)  are  not  accurate  near 

p /p  = 1 .0. 

^en  eb 
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APPENDIX  D.  PROGRAM  DESCRIPTION 


l 

1 

[ A brief  description  of  the  computer  program  is  presented  here  to 

I supplement  the  comments  in  the  code.  Where  appropriate,  the  logic  of 

j each  subroutine  is  discussed.  A list  of  the  programs,  subroutines,  and 

I functions  subprograms  is  presented  in  Table  D-1  . A diagram  of  the  control 

[ logic  followed  by  the  program  is  shown  in  Table  D-2. 

i 

The  computer  program  consists  of  approximately  4000  cards.  It  re- 
quires 30  second  of  central  processor  time  for  compilation  on  the  CDC 
6500  computer.  Nozzle  flow  field  analyses  typically  require  8 to  10 
seconds  of  computation  time,  although  that  time  is  reduced  to  4 to  6 
seconds  when  the  secondary-start  line  option  is  utilized  effectively. 
Boattail  flow  field  analyses  are  substantially  faster  than  nozzle  flow 
field  analyses,  requiring  only  4 to  6 seconds  of  computation  time. 


i 


TABLE  D.l.  LIST  OF  PRORRAI1S,  SUBROUTINES,  AND  FUNCTION  SUBPROGRAMS 


i 


MAIN 

THERMO 

TITLE 

EM 

N0ZZL 

POINT 

BOTTTAL 

INWALL 

NZBSBT 

F 

dPTN 

DRWALL 

0PTB 

WALLPT 

0PTNB 

AXIS 

C0NB0T 

BASE 

PUSH 

TRAN2D 

INPUTT 

TRANl  D 

DATAN 

TRAND4 

DA  TAB 

TRANl  2 

IVLINN 

TRANl  4 

IVLINB 

TRAN24 

XYLINE 

TRAN3D 

STARTN 

TRANSL 

STARTB 

AXITER 

L0GICN 

HILCLM 

L0GICB 

NEWT0N 

WRITE 

GAUSS 

THRST 
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NOZZL 

BOTTAL 

NZBSBT 

OPTN 

OPTB 

OPTNB 

INPUTT 

INPUTT 

INPUTT 

INPUTT 

INPUTT 

INPUTT 

DATAN 

DATAB 

DATAN 

DATAN 

DATAB 

DATAN 

IVLINN 

IVLINB 

IVLINN 

IVLINN 

IVLINB 

IVLINN 

STARTN 

STARTB 

STARTN 

STARTN 

STARTN 

LOGICN 

LOGICB 

L0GICN 

L0GICN 

L0GICN 

DATAB  DATAB 

IVLINB  IVLINB 


STARTB 

L0GICB 

BASE 


MAIN.  Program  MAIN  reads  the  first  data  card,  which  specifies 
whether  an  analysis  or  an  optimization  is  to  he  performed,  and  the 
type  of  geometry.  The  parameter  I0PT  specifies  either  an  analysis 
(I0PT  = 1)  or  an  optimization  ( I0PT  = 2).  ITYPL  specifies  the  type  of 
geometry:  a complete  nozzl e-base-boa ttai  1 system  (ITYPE  = l),a  nozzle  only 
(ITYPE  = 2),  or  a boattail  only  (ITYPE  = 3).  The  first  card  also  contains 
an  alphanumeric  job  title  for  the  run.  When  successive  runs  are  being 
made,  this  data  card  must  be  specified  for  each  run.  Program  MAIN  then 
calls  the  appropriate  analysis  or  optimization  master  control  subroutine. 


TITLE . Subroutine  TITLE  writes  out  the  standard  comments  and  job 
title  on  the  first  page  of  the  output. 


N0ZZL . Subroutine  N0ZZL  contains  the  overall  logic  for  the  nozzle 
flow  fiel'd  analysis.  All  of  the  input  data  for  the  nozzle  are  read  in 
by  calling  subroutine  INPUTT.  Then  subroutine  IVLINN  is  called  to  gener- 
ate the  ini tial -val ue  line,  and  subroutine  STARTN  is  called  to  determine 
the  wall  contour  parameters.  Subroutine  L0GICN,  which  contains  the  logic 
for  determining  the  nozzle  flow  field  by  the  method  of  characteristics, 
is  then  called. 


B0TTAL.  Subroutine  B0TTAL  contains  the  overall  logic  for  the  boat- 
tail  flow  field  analysis.  All  of  the  input  data  for  the  boattail  are 
read  in  by  calling  subroutine  INPUTT.  The  subroutine  IVLINB  is  called  to 
generate  the  initial-value  line,  and  subroutine  STARTB  is  called  to 
determine  the  wall  contour  parameters.  Subroutine  L0GICB,  which  contains 
the  logic  for  determining  the  boattail  flow  field  by  the  method  of 
characteristics,  is  then  called. 


NZBSBT.  Subroutine  NZBSBT  contains  the  overall  logic  for  the  analysis 
of  the  nozzle-base-boa ttail  combination.  There  is  a section  for  the 
nozzle  flow  field  analysis,  a section  for  the  boattail  flow  field  analysis, 
and  a section  for  the  base  region  analysis.  The  analyses  are  performed  in 
the  order  stated.  The  same  logic  pattern  is  followed  for  the  flow  field 
analyses  for  both  the  nozzle  (see  subroutine  N0ZZL)  and  the  boattail  (see 
subroutine  B0TTAL).  The  logic  pattern  reads  in  all  of  the  input  data, 
generates  an  ini tial -val ue  line,  calculates  wall  contour  parameters,  and 
employs  the  method  of  characteristics  to  analyze  the  flow  field.  After 
the  flow  fields  fof  both  the  nozzle  and  boattail  have  been  computed,  sub- 
routine BASE  calculates  the  base  pressure  and  the  base  thrust. 


0PTN.  Subroutine  0PTN  is  the  master  logic  program  for'  nozzle  contour 
optimization.  All  of  the  input  data  are  read  in  and  the  ini tial -val ue 
line  is  generated.  Then  L0GICN  is  called  to  compute  the  secondary-start 
line,  which  is  a right-running  Mach  line  attached  to  the  minimum  allowed 
throat  angle  [ANMIN(l)]  and  extending  either  to  the  axis  of  symmetry  or  to 
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the  exit  plane  of  the  nozzle.  L0GICN  then  restarts  from  this  line,  thus 
increasing  the  efficiency  of  multiple  nozzle  flov;  field  calculations  in  a 
nozzle  or  a nozzle-base-boa ttail  optimization  (see  discussion  in  Section 
II).  The  parameter  IMETH  specifies  vdiich  optimization  algorithm  is  to 
be  employed  (IMETH  = 1,  Axial  Iteration;  IMETH  = 2,  Method  of  Steepest 
Descent;  IMETH  = 3,  Newton's  Method).  The  appropriate  subroutine  is 
then  called  to  perform  the  nozzle  contour  optimization. 


0PTB . Subroutine  0PTB  is  the  master  logic  program  for  boattail 
contour  optimization.  All  of  the  input  data  are  read  in  and  the  initial- 
value  line  is  generated.  The  optimization  method  specified  by  the  para- 
meter IMETH  (see  subroutine  0PTN)  is  called  to  perform  the  boattail 
contour  optimization. 


0PTNB . Subroutine  0PTNB  is  the  master  logic  program  for  a complete 
nozzle-base-boattail  contour  optimization.  All  of  the  input  data  are 
read  in  and  the  initial-value  line  for  the  nozzle  is  generated.  L0GICN 
is  called  to  compute  the  nozzle  secondary-start  line  (see  discussion  in 
Section  II).  The  optimization  method  to  be  employed  is  specified  by  the 
parameter  IMETH  (see  subroutine  0PTN) . The  appropriate  subroutine  is 
called  to  perform  the  combined  nozzle-base-boattail  optimization. 


C0NB0T.  Subroutine  C0NB0T  calculates  the  attachment  angle  corres- 
ponding to  a conical  boattail,  given  the  boattail  exit  radius  and  initial 
expansion  contour.  This  subroutine  employs  a Newton's  method  zero  finding 
algorithm  to  iteratively  solve  for  the  attachment  angle. 


PUSH.  Function  PUSH  calls  the  appropriate  subroutines  to  compute  the 
thrust  for  the  specified  geometry  using  the  previously  specified  start 
lines  for  the  nozzle  and/or  the  boattail.  The  wall  geometry  for  a nozzle- 
base-boattail  configuration  is  specified  by  the  following  three  indepen- 
dent parameters:  (1)  the  nozzle  throat  attachment  angle  Ogp,  (2)  the 

nozzle  exit  radius  yen>  and  (3)  the  boattail  exit  radius  yeb  (see  Figures 
E-1  and  E-2).  Subroutine  PUSH  sets  the  variable  PUSH  equal  to  the  value 
of  the  thrust  for  the  specified  geometry,  and  returns  to  the  calling 
program.  This  function  subprogram  serves  as  the  objective  function  for 
the  optimization  algorithms. 


INPUTT.  Subroutine  INPUTT  reads  in  all  of  the  input  data  through 
NAMELIST  INF0.  If  an  input  variable  is  not  specified  on  the  data  cards, 
then  that  variable  is  set  equal  to  its  default  value  given  in  subroutine 
INPUTT.  A discussion  of  the  input  is  given  in  Appendix  E. 


DATAN . Subroutine  DATAN  defines  all  of  the  variables  employed  in 
the  method  of  characteristics  calculations  to  be  those  corresponding  to 
the  nozzle.  All  of  the  nozzle  input  variables  are  printed  out  at  the  be- 
ginning of  the  analysis. 
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DATAB.  Subroutine  QATAR  defines  all  of  the  variables  employed  in 
the  method  of  charac teristics  calculations  to  be  those  corresponding 
to  the  boattail  . All  of  the  boattail  input  variables  are  printed  out 
at  the  beginning  of  the  analysis. 


IVLINN.  Subroutine  IVLINN  generates  an  ini tial -val ue  line  for  the 
nozzle  flow  field  analysis.  When  the  input  parameter  IVLN  = 1,  the 
points  on  the  initial-value  line  are  read  in.  Each  data  card  must  contain 
the  location,  velocity  magnitude,  pressure,  flow  angle,  and  density  at  a 
point  on  the  i ni tial -va 1 ue  line.  The  number  of  data  cards  must  be  equal 
to  the  parameter  NPN.  If  the  parameter  IVLN  = 0,  Kliegel's  initial-value 
line  is  generated.  This  line  begins  on  the  wall  at  the  minimum  cross 
section  oi  the  throat  and  ends  downstream  a distance  EPSL0N  along  the 
axis.  EPSL0N  is  computed  in  subroutine  IVLINN  and  is  a function  of  the 
specific  heat  ratio,  throat  radius,  and  throat  upstream  radius  of  curva- 
ture. All  of  the  properties  for  each  point  on  the  ini tial -val ue  line 
are  stored  in  the  data  storage  array  SL(K,I)  and  SLN(K,I). 


IVLINB . Subroutine  IVLINB  determines  the  i ni tial -va lue  line  for 
the  boattail.  If  the  input  parameter  IVLB  = 1,  the  points  on  the  initial- 
value  line  are  read  in.  Each  data  card  must  contain  the  location, 
velocity  magnitude,  pressure,  flow  angle,  and  density  at  a point  on  the 
initial-value  line.  The  number  of  data  cards  must  be  equal  to  the  para- 
meter NPB.  If  the  parameter  IVLB  = 0,  a uniform  flow  i ni tial -val ue  line 
is  generated.  This  uniform  line  is  computed  from  the  free-stream  Mach 
number  and  the  stagnation  properties  which  are  specified  in  the  input 
data.  All  of  the  properties  for  each  point  on  the  initial-value  line  are 
stored  in  the  data  storage  arrays  SL(K,I)  and  SLB(K,I). 

Each  point  on  the  i ni tial -val ue  line  is  the  beginning  of  a new  right- 
running Mach  line.  The  ini tial -val ue  line  must  be  specified  long  enough 
to  insure  that  the  last  point  on  the  initial-value  line  will  generate  a 
right-running  Mach  line  beyond  the  exit  lip  point  of  the  boattail.  For 
the  uniform  flow  ini tial -val ue  line,  the  x-coordinate  spacing  of  the 
points  is  determined  by  the  input  parameter  SPACE.  If  SPACE  = 1.0,  the 
x-coordinate  of  the  last  point  on  the  ini tial -value  line  will  be  equal 
to  the  boattail  length  Xgp,  and  the  last  right-running  Mach  line  will 
always  go  beyond  the  wall  exit  lip  point  of  the  boattail.  To  keep  from 
generating  unnecessary  points  on  the  ini tial -val ue  line,  SPACE  should  be 
specified  between  0.3  and  O.C,  depending  on  the  freestream  Mach  number. 


XYL INE . Subroutine  XYLINE  generates  geometrically  progressively 
spaced  y-coordina tes  and  the  corresponding  x-coordinates  for  points 
along  the  nozzle  initial-value  line.  These  points  may  be  closely  spaced 
next  to  the  nozzle  wall  and  spread  apart  as  they  approach  the  nozzle  axis. 
The  geometric  progression  is  determined  by  the  input  parameter  RATI0I. 

For  example,  if  RATI0I  = 5.0,  then  the  y distance  between  the  last  two 
points  on  the  ini tial -val ue  line  at  the  axis  will  be  5.0  times  greater  than 
the  y distance  between  the  first  two  points  at  the  wall.  This  type  of 
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point  spacing  permits  a closer  spacing  of  points  near  the  wall  where 
the  flow  property  gradients  are  larger.  If  RATI01  = 1.0,  the  y- 
coordinates  will  be  uniformly  spaced  along  the  ini tial -value  line. 


S_TA_RTI'l.  Subroutine  STARTN  deteniiines  the  geometric  parameters 
that  specify  the  nozzle  contour.  The  type  of  wall  contour  and  the  para- 
meters that  specify  the  wall  contour  are  printed  out.  A conical  wall, 
a second-order  polynomial  wall,  or  a tabular  wall  contour  may  be 
selected.  The  input  parameter  IWALLN  determines  which  type  of  wall  con- 
tour is  selected.  IWALLN  = 1 specifies  a conical  nozzle,  IWALLN  = 2 
specifies  a second-order  polynomial  nozzle,  and  IWALLN  = 5 specifies  a 
tabular  nozzle.  The  half-angle  of  the  conical  nozzle  is  specified  by 
the  throat  expansion  angle  Oa^i.  For  a second-order  polynomial  wall, 
the  wall  exit  lip  point  is  determined  by  specifying  either  the  exit 
radius  y^p  or  the  exit  angle  Ogn  (see  Figure  E-1).  Subroutine  STARTN 
computes  the  coefficients  for  tne  conical  and  second-order  polynomial 
wall  contours  for  the  following  equation. 

y(x)  = a + bx  + cx^  (D-1 ) 


The  number  of  points  in  a tabular  wall  is  specified  by  the  parameter 
NWALLN,  and  the  arrays  XWN(IOO)  and  YWN(IOO)  contain  the  respective  x- 
and  y-coordinates  of  the  NWALLN  points.  The  conical  wall  option  and 
the  tabular  wall  option  may  be  employed  only  with  the  flow  field 
analysis  option.  For  the  optimization  option,  the  second-order  poly- 
nomial wall  contour  option  must  be  specified. 


STARTS.  Subroutine  STARTS  determines  the  geometric  parameters 
that  specify  the  boattail  contour.  The  type  of  wall  contour  and  the 
corressonding  parameters  are  printed  out.  Four  types  of  boattail 
contours  may  be  selected:  conical,  cylindrical,  second-order  poly- 

nomial, or  tabular  contour.  The  parameter  IWALLB  determines  which  type 
of  wall  contour  is  selected.  IWALLS  = 1 specifies  a conical  boattail, 
IWALLS  = 2 specifies  a second-order  polynomial  boattail,  IWALLS  = 3 
specifies  a cylindrical  boattail,  and  IWALLB  = 5 specifies  a tabular 
boattail.  The  half-angle  of  the  conical  boattail  is  specified  by  the 
initial  expansion  angle  Sab-  a second-order  polynomial  wall,  the 
wall  exit  lip  point  is  determined  by  specifying  either  the  exit  radius 
yeb  or  the  exit  angle  Subroutine  STARTS  computes  the  coefficients 

for  the  boattail  wall  contours  employing  the  following  equation. 

y(x)  = a + bx  + cx^  (D-2) 

The  number  of  points  in  a tabular  wall  is  specified  by  NWALLB,  and  the 
arrays  XWB(IOO)  and  YWB(IOO)  contain  the  respective  x-  and  y-coordinates 
of  the  NWALLB  points.  The  tabular  wall  option  may  be  employed  only 
with  the  flow  field  analysis  option.  For  the  optimization  option,  the 
conical  wall  contour  option  must  be  specified  for  a complete  nozzle- 
base-boattai 1 configuration.  For  the  optimization  of  a boattail  contour 
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by  itself,  the  second-order  quc-viratic  v.'ill  contour  option  is  employed.’ 
When  a complete  nozcle-base-boatta  i 1 anal  is  is  beintj  conducted,  the 
base  pressure  model  is  restricted  to  i-.oattai's  which  have  the  same  x- 
coordinate  for  the  wall  exit  lip  point  and  the  nozzle  exit  lip  point. 
However,  if  anotner  base  pressure  model  is  inserted  whicli  includes  the 
effects  of  different  x-cooroi nates  for  the  exit  lip  points  of  the 
nozzle  and  boattuil,  this  will  in  no  way  affe'  t any  other  part  of  the 
analysis . 


L0GICN.  Subroutine  L0bUJN  contains  the  logic  which  generates  the 
solution  network  for  the  nozzle  flow  field.  Left-running  Mach  lines 
are  followed  to  produce  the  characteristic  mesh  which  constitutes  the 
solution  network.  These  left-running  Mach  lines  originate  from  either 
the  ini tial-val ue  line  or  the  nozzle  axis.  L0GICN  determines  which  of 
these  starting  points  is  to  be  used  for  each  left-running  Mach  line. 

In  addition  to  determining  the  starting  point,  L0GICN  monitors  how 
many  interior  points  are  to  be  calculated,  whether  an  indirect  or  di- 
rect wall  point  is  to  be  computed,  and  when  to  stop  the  calculations. 

The  following  paragrapfis  give  a more  detailed  description  of  the  meqh- 
anics  of  subroutine  L0GICN. 

A thorough  knowledge  of  the  1-J  characteristics  coordinate  system 
shown  in  Fig.  D-l  is  required  to  understand  the  overall  solution  logic. 
Each  I-J  point  represents  a calculated  point  in  the  characteristic- net- 
work. The  I-coordinate  lines  represent  right-running  Mach  lines,  and 
the  J-coordinate  lines  represent  left-running  Mach  lines.  The  point 
(MPN,1)  is  the  first  point  of  the  solution  network. 

After  initiating  various  control  variables,  the  first  point  in  the 
characteristic  network  is  obtained.  As  shown  in  Figure  D-l,  this 
point  (NPM,1)  is  located  at  the  wall  on  the  i ni t i al -val ue  line. 

When  the  main  loop  is  entered,  a decision  must  be  made  immediately 
as  to  which  type  of  starting  point  is  to  be  used.  The  first  NPN  - 1 
times  through  the  loop  an  ini tial -val ue  line  point  is  selected.  There- 
after, an  axis  point  is  selected.  This  axis  point  is  located  at  the 
intersection  of  the  nozzle  axis  and  the  right-running  Mach  line  from 
the  second  point  on  the  previous  left- running  Mach  line.  After  ob- 
taining the  starting  point  for  the  next  left-running  f^ch  line,  interior 
points  are  calculated  along  the  line,  and  the  variable  I is  incremented 
by  one  after  each  point  is  calculated.  Tl"'-  procedure  continues  until 
the  I-coordinate  of  the  interior  point  is  equal  to  ILIMIT,  when  a wall 
point  must  be  computed.  After  the  wall  point  is  coiiiputed,  the  nozzle 
thrust  is  computed  by  subroutine  THRST. 

Tvio  different  methods  are  used  to  compute  a wall  point.  Along  the 
circular  arc  initial  expansion  region  of  the  throat,  where  property 
gradients  are  high,  points  are  specified  along  the  wall  through  which 
the  left-running  Mach  lines  mu'  t pa^s.  Beyond  the  initial  expansion 
region  of  the  circular  arc  throat,  tne  lett-running  Mach  lines  are  ex- 
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Figure  D-1 . Characteristic  Coordinate  System  for  the 
Nozzle  Flow  Field 
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tettded  directly  to  ir.terbect  the  woll.  lln  latter  technique  does  not 
require  any  special  nodi  ti  ea  t;  'H  of  ttu;  loqi..  Put  tne  former  technique 
does.  To  handle  that  situation,  the  S('ecial  ■'uiiiPerinij  system  shown  in 
Fig.  D-2  has  been  devised.  .■  -pecifieo  point,  along  the  wall  are 
spaced  geometrically  to  increa.e  the  accuracy  of  the  -.olution.  Ttie 
parameter  ItATiyiW  specifies  tne  ratio  of  the  ainjuldr  iiicrement  spanned 
by  the  last  two  points  on  the  initial  exfiansion  lontour  adjacent  to 
the  throat  attaciiment  angle  to  the  angular  increment  spanned  by  the 
first  two  points  at  the  throat.  Knowimi  the  throat  attachment  RATI0W, 
the  location  of  the  throat,  and  the  nozzle  downstream  radius  of 
curvature,  the  wall  point  locations  may  be  computed.  The  interpolated 
points  labeled  by  I's  on  Figure  D-r  re  not  stored.  The  calculated  wall 
point  (I  ^ 1,  J - 1]  ill  .'trated  in  ! inure  l'-2(a}  is  stored  in  the 
previous  1 ef t - rynni ng  Mach  li’io  array  at  the  (I  + 1)  location.  After 
storing  and  priming  the  wall  poiiit,  th.f‘  interior  point  (1  '■  ] , j) 
illustrated  in  Mgnre  U-2(h)  is  calculated.  This  entire  procedure, 
which  is  referred  co  as  the  inverse  wall  point  scheme,  is  repeated  until 
a specified  wall  point  requires  the  ;<-value  of  the  interpolated  point 
1 to  exceed  that  ol  ttie  last  inferior  point.  When  this  occurs,  the 
parameter  ITER  is  set  equal  to  9'1  in  subroutine  INWAll.,  and  the  next 
left-running  Mach  line  is  scarpno  from  either  the  initial-value  line  or 
the  axis.  Tne  lasi.  invei  .e  wall  point  is  computed  at  the  maximum  throat 
angle  Gap. 

The  main  loop  continues  to  generate  left-running  Mach  lines  by  the 
direct  wall  point  nechod  until  the  x-coordinate  of  the  wall  point  ex- 
ceeds tlie  specified  nozzle  length  x^p.  Then  the  inverse  vwll  point 
routine  is  called,  using  the  nozzle  exit  lip  point  as  the  specified  wall 
point.  This  completes  the  nozzle  flow  field  analysis. 

After  the  calculation  of  each  interior  point,  a check  is  made  to 
see  if  the  Mach  lines  crossed.  If  the  right-running  Mach  lines  crossed, 
then  the  point  just  calculated  is  dropped  from  the  flow  field  and  a 
message  is  printed  indicating  this  has  beer,  done.  The  properties  at  the 
next  interior  point  are  then  based  on  the  next  point  along  the  previous 
left-running  Mach  line.  If  the  left-running  Mach  lines  cross,  the  cal- 
culation of  points  along  the  new  left-running  Mach  line  is  terminated. 

The  remaining  point;,  up  to  the  wall  are  defined  to  be  those  of  the 
previous  left-running  Mach  line  which  was  just  crossed.  A message  is 
printed  indicating  tiiat  the  left-running  Mach  lines  crossed. 

Two  inipnrtint  storage  procedures  are  carrier!  our  in  subroutine 
L0GICN.  A secondary-start  line,  which  is  eiirmoyed  during  the  optimiza- 
tion procedure,  is  stored  when  the  paramerei  I'.SL  is  non-zero  (see 
Section  II).  The  secondary-start  line  attaches  to  the  nozzle  vjal1  at 
the  minimum  allov/ed  throat  attach'iient  angle.  It  is  stored  in  the  same 
location  as  the  initial-value  line.  Subroutine  10GICN  restarts  the 
analysis  from  tiiis  secondary-s r.a rt  line.  The  value  for  ISSL  = NPN  + 
riPW  when  the  secondary-st.'irt  line  is  being  stored.  When  the  complete 
seconda ry-star*  line  I, a;  been  detepniivd,  the  value  of  ISSL  is  0.  Note 
that  subsequent  restarts  mir  ‘ o nreiter  throat  attachment  angles  than 
the  min-imum  angle  v.'hicii  deterniir.'-'d  t^‘i  secondary-start  line.  If  the 
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Figure  D-2.  Characteristic  Coordinates  for  an  Inverse  Wall 
Point  and  an  Interior  Point 
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dnalyiis  option  is  jtiliz-.  O,  IS'-t  i-  .oi  c to  a laryo  number  and  a 
secondary-star;,  line  is  tan;  :.  red.  Tne  i it  properties  at  the  final 
wall  point  of  the  flow  fi' ’ -e  stered  ''n  the  array  Dri'lAL. 


L0GICB.  sLibt  cutii.L-  i.dulii;  cuniJiris  th'.-  k'/i.  that  generates  the 
solution  network  for  ttie  Li-j.ittail  flow  field.  Pi ijht-runnincj  Mach 
lines  are  followed  f--oi:;  the  ini  t i a 1 -va  1 ue  line  to  produce  the  character- 
istic meslu  bubn.ul  ii\e  determines  how  man;,'  interior  points  are 

to  be  calculated,  when  a ’wall  ,.<oh'il  is  to  be  coiiiputed,  and  when  to  stop 
the  calculations.  Hie  fonowir!i  paragraph.,  liive  a detailed  description 
of  the  mechanic,;,  :f  ;>ut''Ojtine  l,'blf['. 

As  in  the  ne,cle  i b ..  ;u'id  ..inul/sis,  an  i-J  cnar.icter i Stic  coordi- 
nate systair.  is  e!.plo,eci  to  perfoma  tiie  crackiiu;  and  storage  of  tfie  com- 
puted points  in  the  boattail  t’ow  field.  Tne  coordinate  system  is 

shown  in  f^ig.  D-,>.  The  i -coordinate  lines  represent  right-running 
Mach  lines,  and  tluj  j-eoordi  nato  li'Mis  represen*  left-running  Mach  lines. 
Ihe  point  (1,1)  is  the  fir->i.  print  >.:f  t!ie  solution  network. 

Tne  main  los , ot  rl  rcu*  ire  lidCidF  •*„  exceptionally  straightforward 
since  each  ri ;,h*. -r.;nni >v;  I', a..;*,  line-  i;'.ust  originate  from  the  initial-value 
line  and  tenninate  at  tse  w-;ll.  Thus,  rhe  first  step  in  the  loop  is  to 
obtain  a point  ■•'rom  the  in  i tiai  - v : 1 ue  line.  Inserior  point  calculations 
are  then  performed  until  tiih  .‘-coordinate  equals  JLIMIT.  Then  a direct 
v/a1l  point  is  deten-ined,  and  th-;  thrust  associated  with  that  point  is 
computed.  The  par.j.neter  JL ''.1'^  i set  equal  to  d a-  1 , and  the  entire 
loop  is  repeated,  'when  ■ ‘le  — coordin  ;le  ot  the  wal’  point  exceeds  the 
boattail  lemith  xop  the  invc'se  wall  point  routine  is  called  to  deter- 
mine the  flow  propertie;.  at  the  boattail  exit  lip  point.  The  thrust  is 
then  computed.  The  :inal  pome  on  the  last  Mach  line  is  stored  in  the 
array  DFIhAL..  This  cosiuleses  she  ilow  fie^d  analysis  for  the  boattail. 

If  the  i ni tiai -val ue  line  for  the  boattail  flow  field  analysis  has 
been  specified  too  s'  irt,  tne  ..uvil/sis  wik  be  contiued  until  the  last 
poinc  on  the  i n i ti  .-1  - ,a  i ut  lio.e  ‘.j:  beeii  used,  at  which  time  the  last 
wall  point  still  loes  i:ot  eAC-md  the  S’pecified  boattail  length  Xen- 
When  this  occurs,  an  erro-'  ; ;;  ‘.-rtne  will  f.e  printed  and  the  program  stops. 


WJkTT  Subroutine  'St'*  '*ut  tiwi:  resuits  of  the  character- 

istic alculusions  t'  ■ ; ; • ’ r.dar'-'.  ta.  ‘ ‘in'.-.  !t  also  detemines  the 
prooer  neadk';:>  an;  ;. .. 'I'iu  .=  -.ontr;ds  ner.'  -'Miry  for  the  various  print 
options. 


THRST.  Subr;,./  TiiK..’  .-Iciik’tes  the  thrusi.  at  a wall  point  for 
the  noczl-..-.  for  so,  ,,a .'-i  . ;.;;e  one-dr-ie:'  'oii  il  flow  thrust  is  also 

computed.  This  onc-dimensinn.k  thrust  aopears  on  the  last  line  of 
output  for  each  I c-f t-runniii.s;  M-ich  i me  in  tne  nozzle.  Subroutine  THRST 
also  computes  th-  ■‘"rmris--,  ! 'le  *h'mt  and  th‘>  mass  flow  rare  between 

an  interior  poms  and  Inr  i.-i 


INITIAL- VALUE 
LINE  ^ 


Figure  D-3.  Characteristic  Coordinate  System  for  the 
Boattail  Flov;  Field 


T_HERM0.  Subroutine  THLRfW  calculdtes  •.i:e  thermodynamic  properties 
of  the  gas  which  are  required  by  tt'O  charactei  istic  and  compatibility 
equations.  The  Mach  number.  f>’,  temperature,  t,  and  the  speed  of  sound, 
a,  are  computed  as  a function  of  the  density,  , . pressure  p,  and  velocity 
V.  The  subroutine  built  into  the  program  is  based  on  the  properties  of 
a thermally  and  calorically  perfect  gas.  Other  gas  thermodynamic  nx)dels 
may  be  considered  by  appropriate  modifications  to  subroutine  THERM0. 


EM.  Function  EM  determines  the  Mach  number  as  a function  of  the 
area  ratio  for  the  one-dimensional  flow  of  a perfect  gas.  That  informa- 
tion is  required  in  subroutine  THRST  for  ttie  one-dimensional  thrust 
comparison.  !f  a gas  thermodynami c model  other  than  a perfect  gas  is 
specified  in  subroutine  TflERM0,  function  EM  must  be  modified  accordingly. 


W^NT.  Subroutine  R0IRT  evaluates  the  solution  at  an  interior 
point  in  the  flow  field  using  the  method  of  characteristics  for  rota- 
tional flow.  The  solution  process  employs  trie  modified  Euler  predictor- 
corrector  technique  for  solving  the  flow  field  characteristic  and 
compatibility  equations.  An  average  property  scheme  is  used  to  compute 
the  coefficient  of  each  differential.  The  point  labeling  scheme  for 
this  subroutine  is  shown  in  Fig.  D-4.  The  location  and  flow  properties 
are  known  at  points  1 and  2.  Point  4 is  located  at  the  intersection  of 
the  right-running  and  left-running  Mach  lines.  The  location  of  point  3 
is  obtained  by  extending  a streamline  back  from  point  4 until  it  inter- 
sects the  line  between  points  1 and  2.  Flow  properties  at  point  3 are 
obtained  by  linear  interpolation  between  points  1 and  2.  The  flow 
properties  p,  . , V,  and  at  point  4 are  then  obtained  by  solving  the 
characteristic  and  compatibility  equations.  The  equations  programmed 
are  as  follows.  Along  the  right-running  Mach  line. 




^ ^ yT-r-cos^”-)  ^ 0 


Along  the  left-running  Mach  line. 


-dx  = ^ 
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dp  - d - r ; qy  = Q 


(D-3) 

(D-4) 


(D-5) 

(D-6) 
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STREAMLINE 


RIGHT-RUNNING 
MACH  LINE 


Figure  D-4.  Point  Labeling  Scheme  for  an  Interior  Point 


Along  the  streamline. 


dy^  _ y^ 
dx  u 

(D-7) 

pVdV  + dp  = 0 

(D-8) 

2 

dp  - a dp  = 0 

(D-9) 

The  modified  Euler  predictor-corrector 
is  used  in  the  numerical  analysis.  The 

technique  with  average  properties 
remaining  flow  properties  are  ob- 

tained  by  calling  subroutine  THERM0. 


INWALL.  Subroutine  INWALL  evaluates  the  solution  at  a specified 
wall  point  in  the  circular  arc  throat  region  using  an  inverse  method  of 
characteristics  for  rotational  flow.  The  location  and  flow  properties 
are  known  at  points  2,  3,  and  5 in  Fig.  D-5.  The  solution  is  sought  at 
point  4 where  the  location  (x4,y4)  and  slope  64  are  already  known.  A 
left-running  Mach  line  is  extended  rearward  from  point  4 to  point  1, 
where  the  flow  properties  are  obtained  by  linear  interpolation  between 
points  2 and  3.  The  compatibility  equations  are  then  solved  for  the 
flow  properties  p,  p,  V,  and  0 at  point  4.  The  equations  programmed 
are  as  follows.  Along  the  left-running  Mach  line. 


^ = tan  (0  + a)  (D-10) 


I 


1 


dp  - do  + 


6 [ 


V 

yM\T  sin  (0  + a) 


dy  = 0 


(D-11) 


Along  the  nozzle  wall  , 


dx 

V 

u 

(D-12) 

pVdV 

+ dp  = 0 

(D-13) 

dp  - 

a^do  = 0 

(D-14) 
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CIRCULAR  ARC  CONTOUR  OF 


Figure  D-5.  Point  Labeling  Scheme  for  an  Inverse  Wall  Point 


The  niodified  Euler  oredi ctor-corrector  technique  with  average  properties 
is  used  in  the  numerical  analysis.  The  remaining  flow  properties  are 
obtained  by  calling  subroutine  THERM0. 


F_.  Function  F calculates  the  nozzle  wall  point  spacing  required 
for  the  portion  or  the  circular  arc  throat  between  the  secondary-start 
line  attachment  point  and  the  circular  arc  throat  attachment  point.  The 
objective  is  to  maintain  consistent  spacing  for  those  points  upstream 
of  the  secondary-start  line  attachment  point  and  those  downstream  of 
that  point. 

First,  function  F computes  the  number  of  wall  points  from  the  mini- 
mum throat  attachment  point  to  the  specified  throat  attachment  point 
using  the  original  wall  point  spacing  between  the  nozzle  throat  and  the 
minimum  throat  attachment  point.  Then  an  iteration  is  performed  to 
find  the  exact  spacing  which  allows  the  last  point  to  lie  at  the  speci- 
fied throat  attachment  point.  Thus,  the  inverse  wall  points  located 
downstream  of  the  secondary  start  line  point  are  spaced  consistently 
with  those  upstream  of  that  point,  and  the  last  inverse  wall  point  is 
the  specified  throat  attachment  point. 


b 


DRWALL.  Subroutine  DRWALL  determines  the  solution  at  a direct 
wall  point  using  the  method  of  characteristics  for  rotational  flow.  The 
location  and  flow  properties  are  known  at  points  2 and  3 in  Fig.  D-6. 
Since  the  wall  contour  y(x)  and  the  wall  slope  dy(x)/dx  = tanG  are  known, 
it  is  only  necessary  to  locate  point  4 on  the  wall  and  to  calculate  the 
flow  properties  p,  p,  and  V at  that  point.  Left-running  characteristics 
and  compatibility  equations  are  used  for  the  nozzle  wall,  and  right- 
running characteristic  and  compatibility  equations  are  used  for  the 
boattail  wall.  The  equations  programmed  are  as  follows.  Along  the 
Mach  lines, 


= 

dx 

tan  (0  ± a) 

(D-15) 

pv^ 

dp. 

T dO^  + 

yMV  sin  (o  ± ct) 

] dy. 

= 0 

■i 

(D-16)  j 

The  upper  si 

gns 

are  used 

i for  the  nozzle  ana 

lysis. 

and  the 

■j 

lower  signs  ] 

are  used  for 

the 

boatta i 

1 analysis.  Along 

the  wa 

11s, 

dy.  = 

dx 

V 

u 

(D-17)  1 

j 

oVdV 

+ dp  = 0 

(D-18)  1 

dp  - 

a^dp  - 0 

(D-19)  j 

r 


Figure  D-6.  Point  Labeling  Scheme  for  a Direct  Wall  Point 


The  ntodified  Tuler  pre'ii ctor-corrector  tedinique  vjith  average  properties 
is  employed  in  the  numerical  analysis.  T -e  remaining  flow  properties 
are  obtained  by  calling  subroutine  THTRI'10. 


WALLPT.  Subroutine  UALLPl  locates  the  intersection  of  a Mach  line 
with  the  wall  contour.  The  x-coordinate , y-coordinate,  and  wall  slope 
are  computed  for  the  solution  point  and  returned  to  subroutine  DRWALL. 
If  a tabular  wall  has  been  specified,  then  conical  wall  segments  are 
fit  between  each  two  data  points  on  the  wall. 


AXIS.  Subroutine  AXIS  determines  the  solution  at  an  axis  point 
using  the  method  of  characteristics  for  rotational  flow.  The  location 
and  flow  properties  are  known  at  points  1 and  3 in  Fig.  D-7.  Since  the 
flow  is  axisymmetric , the  flow  angle  0 and  y-coordinate  at  the  axis 
must  be  zero.  The  compatibility  equations  along  the  right-running  Mach 
line  and  the  axis  are  solved  to  obtain  the  flow  properties,  p,  p,  and 
V.  The  equations  programmed  are  as  follows.  Along  the  Mach  line, 

= tan  (O  - x)  (D-20) 

dp  + dG  + 6 [ ^mTcosTg^-'^O-^  ° 

Along  the  axis , 

pVdV  + dp  = 0 (D-22) 

dp  - a^dp  = 0 (D-23) 

The  modified  Euler  predictor-corrector  technique,  based  on  empolying 
average  properties  for  the  corrector,  is  used  in  the  numerical  analysis. 
The  remaining  flow  properties  are  obtained  by  calling  subroutine  THERM0. 


BASE.  Subroutine  BASE  calculates  the  base  pressure,  the  base 
thrust,  and  the  total  nozzle-base-boattail  thrust.  The  exit  lip  point 
properties  for  the  nozzle  and  the  boattail  arc  known  from  the  nozzle 
and  boattail  flow  field  calculations.  Using  the  known  exit  lip  point 
properties  in  the  oase  pressure  model  derived  from  the  curve  fit  of 
Addy's  (3)  data,  the  base  pressure  and  the  base  thrust  are  calculated. 
The  total  thrust  is  summed  and  printed  out.  Subroutine  BASE  also  notes 
if  the  boattail  exit  radius  is  less  than  the  nozzle  exit  radius,  a 
physical  impossibility.  A message  to  this  effect  is  printed  out,  and 
the  program  continues. 


So 


LEFT-RUNNING 
MACH  LINE 


Figure  D-7.  Point  Labeling  Scheme  for  an  Axis  Point 
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T^_N2_D.  Subroutine  TRAN2D  transfers  the  properties  stored  at 
location  D (I,J,N)  in  the  D array  to  point  2. 


IRANI D.  Subroutine  IRANI D transfers  the  properties  stored  at 
location  D (I,J,N)  in  the  D array  to  point  1. 


IRAND4.  Subroutine  IRAND4  transfers  the  properties  stored  at 
point  4 to  the  D array  at  location  D (I,J,N). 


IRANI 2.  Subroutine  IRAN12  transfers  the  properties  stored  at 
point  2 to  point  1 . 


IRANI 4.  Subroutine  IRAN14  transfers  the  properties  stored  at 
point  4 to  point  1 . 


IRAN24.  Subroutine  IRAN24  transfers  the  properties  stored  at 
point  4 to  point  2. 


IRAN3D.  Subroutine  IRAN  3D  transfers  the  properties  stored  at 
location  D (I,J,N)  in  the  D array  to  point  3. 


IRANSL . Subroutine  IRANSL  transfers  the  properties  stored  at 
point  4 to  the  secondary-start  line  array  SL  (K,JJP). 


AXIIER.  Subroutine  AXIIER  performs  an  axial  iteration  optimiza- 
tion of  the  objective  function,  thrust  (PUSH).  Ihe  subroutine  requires 
an  initial  estimate  (ANGA,  RADB,  ANGC,  RADD)  and  increments  (DA,  DB, 

DC,  DD)  for  each  of  the  independent  variables.  Ihe  variables  ANGC  and 
DC  are  employed  only  when  optimizing  a boattail  contour  by  itself.  It 
then  perturbs  each  of  the  variables  in  turn,  fitting  a quadratic  poly- 
nomial of  the  form  y = a + bx  + cx^  along  the  variable  direction.  Ihe 
algorithm  then  finds  the  maximum  point  of  the  quadratic  polynomial,  cal- 
culates the  objective  function  there,  and  steps  to  the  maximum  of  the 
points  on  this  line.  Ihe  algorithm  then  proceeds  to  the  next  indepden- 
dent  variable.  Note  that  since  all  of  the  variables  are  independent, 
the  search  directions  are  orthogonal.  If  the  algorithm  encounters  a 
prespecified  boundary  in  any  of  the  variables,  this  is  noted  and  the 
variable  is  set  equal  to  the  boundary  value.  If  the  algorithm  encounters 
an  inequality  constraint,  such  as  the  nozzle  and  boattail  exit  lip 
points  coming  closer  than  some  previously  specified  value  Ay^,  the 
variables  become  dependent  (yg[j  = yen  Ayb)  variable  space  is 

diminished  by  one  dimension. 


57 


The  optimization  may  then  proceed  until  convergence,  when  the  function 
changes  less  than  some  specified  relative  tolerance.  This  algorithm 
stops  in  some  (ITER0)  number  of  base  point  moves  with  an  error  message. 


HILCLM.  Subroutine  HILCLM  performs  an  optimization  using  the 
method  of  steepest  descent  (hill  climbing).  Beginning  with  an  initial 
estimate  (ANGA,  RADB,  ANGC,  RADD)  and  increments  (DA,  DB,  DC,  DD) , 
a first-order  approximation  to  the  unit  gradient  is  calculated.  The 
base  point  is  then  successively  moved  along  the  direction  of  the  gradient 
until  the  objective  function  decreases.  Then,  the  algorithm  recalculates 
the  gradient  and  begins  again.  This  continues  until  the  first  step 
along  a gradient  no  longer  increases  the  objective  function.  Then  the 
finite  difference  increments  and  step  length  are  halved  and  the  al- 
gorithm continues.  Convergence  is  attained  when  the  magnitude  of  the 
gradient  becomes  less  than  some  specified  relative  tolerance.  Some 
maximum  number  of  gradient  evaluations  (ITER0)  stops  the  algorithm 
with  an  error  message. 


NEWT0N.  Subroutine  NEWT0N  performs  an  optimization  employing 
Newton's  method.  The  algorithm  requires  a starting  point  (ANGA,  RADB, 
ANGC,  RADD)  and  increments  (DA,  DB,  DC,  DD)  to  form  a matrix  of  second 
derivative  approximations  (Hessian)  and  a column  vector  along  the  nega- 
tive gradient  (see  discussion  in  Section  II).  Subroutine  GAUSS  is  em- 
ployed to  invert  the  Hessian  and  solve  for  the  base  point  move,  both 
in  distance  and  direction.  If  this  move  exceeds  the  specified  optimi- 
zation boundaries  (ANMAX,  ANMIN),  or  the  step  does  not  increase  the 
objective  function,  then  some  relaxation  factor  is  employed  which  forces 
the  step  to  remain  in  bounds.  This  then  becomes  the  new  base  point, 
and  the  process  repeats.  Convergence  is  attained  when  the  function 
changes  less  than  some  relative  tolerance.  If  the  algorithm  does  not 
converge  within  ITER0  base  point  moves,  an  error  message  is  printed  and 
the  program  stops. 


GAUSS.  Subroutine  GAUSS  is  a simultaneous  linear  equation  solver. 

A system  of  N linear  equations,  with  coefficient  matrix  C and  right  hand 
side  vector  B is  solved.  The  solution  is  the  column  vector  V.  The 
algorithm  employs  Gauss  reduction,  requiring  only  N - 1 passes  to 
complete  the  reduction.  Back  substitution  then  fills  the  V vector  be- 
fore it  returns  to  subroutine  NEWT0N,  the  calling  program. 
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APPENDIX  L.  INPUT  PAPAflETERS 


The  first  data  card  in  each  run  specifies  the  type  of  computation. 
The  variable  I0PT  contained  in  column  1 selects  whether  an  analysis 
(I0PT  = 1)  or  an  optimization  (I0PT  = 2)  is  performed.  The  variable 
ITYPE  specified  in  column  2 determines  the  type  of  geometry  being  con- 
sidered; ITYPE  = 1 for  a complete  nozzle-base-boattai 1 assembly, 

ITYPE  = 2 for  a nozzle  only,  and  ITYPE  = 3 for  a boattail  only.  The 
remaining  78  spaces  are  available  for  any  alphanumeric  job  title,  which 
is  printed  at  the  beginning  of  each  run. 

All  other  input  parameters  are  read  in  by  subroutine  INPUTT  through 
a NAMELIST  file  called  INF0.  Each  of  the  parameters  has  a default  value 
defined  within  the  program.  All  of  the  input  parameters  necessary  to 
perform  nozzle,  boattail,  or  nozzl e-base-boattai 1 analyses  and  optimiza- 
tions are  included  in  NAMELIST  INF0. 

The  code  has  the  capability  of  using  tabular  values  for  gas  proper- 
ties, wall  contours,  and  i ni ti al - val ue  lines.  If  tabular  gas  properties 
are  not  input,  the  program  uses  a constant  specific  heat  ratio  and  gas 
constant  (table  of  1 row).  However,  to  use  the  tabular  gas  property 
option,  the  number  of  rows  in  the  gas  table  NGASN  or  NGASB,  respectively, 
(up  to  20)  must  be  specified.  Each  row  includes  a velocity  magnitude, 
a specific  heat  ratio,  and  a gas  constant. 

Normally,  the  tabular  wall  contour  option  is  not  employed  and  the 
wall  contour  is  internally  computed  as  either  a conical  contour  or  a 
quadratic  polynomial  contour.  When  the  tabular  wall  option  is  employed, 
the  number  of  points  on  the  tabular  wall  is  given  by  NWALLN  and  NWALLB 
for  the  nozzle  and  boattail,  respectively.  Each  point  must  have  an  x- 
and  y-coordi Mate , up  to  a maximum  of  100  wall  points. 

When  the  tabular  initial-value  line  option  is  not  utilized,  an 
initial-value  line  is  internally  generated,  using  Kliegel's  transonic 
analysis  for  the  nozzle  and  a uniform  flow  for  the  boattail.  The  pro- 
cedure for  specifying  tabular  i ni ti al -val ue  lines  is  discussed  in  the 
descriptions  of  subroutines  IVLINN  and  IVLINB  in  Appendix  D. 

Tables  E-1  through  E-3  list  the  program  input  variables,  a descrip- 
tion of  each  variable,  the  engineering  symbol  for  each  variable,  and 
the  units  and  default  value  of  each  variable.  When  no  units  are 
specified,  the  variable  is  dimensionless.  Figures  E-1  and  E-2  illustrate 
the  nozzle  and  boattail  geometries. 
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TABLE  E.l.  INPUT  PARAMETERS  FOR  NOZZLE 


Variable  Description 

Units 

Defaul  t 
Val  ue 

NPN 

Number  of  points  on  the  initial -val ue  line. 

11 

IVLN 

Type  of  i ni ti al -val ue  line:  0 Klieqel, 
1 read  in  points. 

0 

RAT  10 1 

Geometric  ratio  of  point  spacing  along  the 
initial-val ue  line. 

5.0 

IVLSN 

Terminate  analysis  after  nozzle  initial-value 
line  calculation;  0 no,  1 yes. 

0 

NPW 

Number  of  points  along  throat  initial  expansion 
contour. 

22 

RATI0W 

Geometric  ratio  of  point  spacing  along  throat 
expansion  contour. 

5.0 

ICN 

Number  of  applications  of  the  corrector. 

2 

T0LN 

Fractional  convergence  tolerance  of  the 
corrector. 

0.001 

KWRITN 

Output  option:  0 every  point,  1 initial-value, 
axis,  and  wall  points,  2 wall  points, 

3 exit  lip  point. 

0 

IDN 

Dump  flag:  0 no,  1 yes. 

0 

YTN 

Throat  radius,  y^. 

in. 

1.0 

THETTN 

Throat  attachment  angle,  G,„. 

an 

deg 

30.0 

RUN 

Nozzle  throat  upstream  radius  of  curvature, 
Ptu- 

in. 

3.0 

RDN 

Nozzle  throat  downstream  radius  of  curvature, 
Ptd- 

in. 

0.5 

XEN 

Nozzle  exit  length,  x^^. 

in. 

10.0 

YEN 

Nozzle  exit  radius,  y^j^. 

in. 

5 ■: 

THETEN 

Nozzle  exit  angle, 

an 

deg 

0.0 

IWALLN 

Wall  contour:  1 cone,  2 second-order  poly- 
nomial, 5 tabular. 

2 
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TABLL  t.l.  CONTItWLD 


Variable  Description 

Units 

uc  1 a u t 

Va  ljie__ 

DTHETN 

Incremental  angle  change  in  Prandtl -Meyer 
expansion. 

deg 

1 .0 

DEL TAN 

1.0  axisyiiuiietric  nozzle,  0.0  planar  nozzle. 

1.0 

P0N 

Stagnation  pressure,  P^. 

psia 

1000.0 

T0N 

Stagnation  temperature,  T^. 

R 

6000.0 

PAMB 

Ambient  pressure, 

psia 

0.0 

NGASN 

Number  of  rows  in  gas  property  table, 
if  1,  then  constant  y and  R. 

1 

VEIN 

Tabular  values  of  velocity  magnitude,  V, 
NGASN  values. 

ft 

sec 

none 

DAMN 

Tabular  values  of  specific  heat  ratio,  y, 
NGASN  values. 

1.2 

RGASN 

Tabular  values  of  gas  constant,  R, 
NGASN  values. 

ft-lbf 
1 bm-R 

60.0 

NWALLN 

Number  of  points  on  tabular  wall. 

none 

XWN 

x-coordinates  of  points  on  tabular  wall, 
NWALLN  values. 

in. 

none 

YWN 

y-coordinates  of  points  on  tabular  wall, 
NWALLN  values. 

in. 

none 

TABLE  E.2.  INPUT  PARAMETERS  FOR  BOATTAIL 


Variable  Description 

Units 

Default 
Val  ue 

NPB 

Number  of  points  on  the  i ni tial -val ue 
line. 

30 

IVLB 

Type  of  i ni tial -val ue  line:  0 uniform, 
1 read  in  points. 

0 

IVLSB 

Terminate  boattail  analysis  after  boat- 
tail  initial-value  line  calculation: 
0 no,  1 yes . 

0 

ICB 

Number  of  applications  of  the  corrector. 

2 

T0LB 

Fractional  convergence  tolerance  of 
the  corrector. 

0.001 

KWRITB 

Output  option:  0 every  point,  1 initial- 
value,  axis  and  wall  points,  2 wall 
points,  3 exit  lip  point. 

0 

IDS 

Dumo  flag:  0 no,  1 yes. 

0 

MS 

Freestream  Mach  number. 

1.4 

SPACE 

Spacing  factor  for  points  along  the 
initial-value  line. 

0.5 

YTB 

Initial  boattail  radius,  y^. 

in. 

8.0 

THETTB 

Boattail  attachment  angle,  8^. 

deg 

2.0 

RDB 

Boattail  downstream  radius  of  curvature, 
f^td- 

in. 

5.0 

XEB 

Boattail  exit  length,  x^^. 

in. 

10.0 

YEB 

Boattail  exit  radius,  y^^^. 

in. 

7.0 

THETEB 

Boattail  exit  angle,  0^^. 

deg 

0.0 

IWALLB 

Wall  contour:  1 cone,  2 second-order 

polynomial,  3 cylindrical,  5 tabular 

2 

DTHETB 

Incremental  angle  change  in  Prandtl -Meyer 
expansion . 

deg 

1.0 

TABLE  E.2.  COmiNUED 


Variable  Description 

Units 

Default 

Value 

DEL TAB 

1.0  axi syinnetri c boattail,  0.0  planar 
boattai 1 . 

1.0 

P0B 

Stagnation  pressure,  P^. 

psia 

12.0 

T0B 

Stagnation  temperature,  1^. 

R 

600.0 

NGASB 

Number  of  rows  in  gas  property  table, 
if  1,  then  constant  y and  R. 

1 

VELB 

Tabular  values  of  velocity  magnitude, 
V,  NGASB  values. 

ft 

sec 

none 

GAMB 

Tabular  values  of  specific  heat  ratio 
, , NGASB  val ues . 

1.4 

RGASB 

Tabular  values  of  gas  constant,  R, 
NGASB  values. 

ft-lbf 

Ibm-R 

53.3 

NWALLB 

Number  of  points  on  tabular  wall. 

none 

XWB 

/-coordinates  of  points  on  tabular  wall, 
NWALLB  values. 

in. 

none 

YWB 

y-coordinates  of  points  on  tabular  wall, 
NWALLB  values. 

in. 

none 

TABLE  E.3.  INPUT  PARAMETERS  FOR  OPTIMIZATION 


Variable 

Description 

Units 

De  faul  t 
Value 

ANSTRT(l) 

Initial  estimate  of  nozzle 

throat  attachment  angle,  0 
^ an 

deg 

30.0 

ANSTRT(2) 

Initial  estimate  of  the  nozzle  exit 
radius, 

in. 

5.0 

ANSTRT(3) 

Initial  estimate  of  the  boattail 
attachment  angle,  Oab*  Only  employed 
when  designing  a boattail  alone. 

deg 

15.0 

ANSTRT(4) 

Initial  estimate  of  the  boattail  exit 
radius,  ygb-  Only  employed  when  de- 
signing a boattail  alone. 

in. 

4.0 

ANMAX(l) 

Maximum  allowed  nozzle  throat  attach- 
ment angle. 

deg 

45.0 

ANMAX(2) 

Maximum  allowed  nozzle  exit  radius. 

in. 

7.0 

ANMAX(3) 

Maximum  allowed  boattail  attachment 
angle. 

deg 

20.0 

ANMAX(4) 

Maximum  allowed  boattail  exit  radius. 

in. 

8.0 

ANMIN(l) 

Minimum  allowed  nozzle  throat  attachment 
angle. 

deg 

20.0 

ANMIN(2) 

Minimum  allowed  nozzle  exit  radius. 

in. 

4.0 

ANMIN(3) 

Minimum  allowed  boattail  attachment 
angle. 

deg 

0.0 

ANMIN(4) 

Minimum  allowed  boattail  exit  radius. 

in. 

4.0 

IMETH 

Optimization  method:  1 axial  iteration, 

2 method  of  steepest  descent,  3 Newton' 
method. 

s 

3 

ITER0 

Maximum  number  of  passes  through  the 
optimization  algorithm. 

10 

T0L0 

Relative  convergence  tolerance  for 
optimization 

0.001 
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TABLE  E.3.  COIiTHJUED 


Variable 

Description 

Units 

Defaul t 
Value 

KWRIT0 

Output  option;  0 every  point,  1 
initial-value,  axis,  and  wall 
points,  2 wall  points,  3 exit  lip 
point. 

3 

I0P 

Specifies  type  of  geometry  for  optimi- 
zation: 1 nozzle,  2 boattail,  3 
nozzle-base-boattai 1 combi  nation. 

1 

N 

Specifies  the  dimension  of  the  optimi- 
zation space  (i.e.,  the  number  of 
unconstrained  parameters). 

2 

DA 

Nozzle  attachment  angle  perturbation. 

deg 

1 .0 

DB 

Nozzle  exit  lip  radius  perturbation. 

in. 

0-2ytn 

DC 

Boattail  attachment  angle  perturba- 
tion. 

deg 

0.5 

DD 

Boattail  exit  lip  radius  perturba- 
tion. 

in. 

0-2ytn 

DBAS 

Minimum  annular  base  width,  Ay^. 

in. 

0.0 
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Specification  of  the  Boattail  Geometry 


APPENDIX  F. 


SAMPLE  CASES 


<^ive  sample  cases  are  presented  in  this  section  to  illustrate  the 
application  of  the  computer  program.  These  cases  exercise  most  of  the 
options  available  for  analysis  or  optimization. 

Each  data  deck  begins  with  a TITLE  CARD.  That  card  identifies  the 
type  of  geometry  (nozzle,  boattail,  or  nozzle-base-boattail  assembly) 
and  whether  an  analysis  or  an  optimization  is  desired.  Column  1 con- 
tains the  parameter  I0PT  which  is  1 for  an  analysis  and  2 for  an 
optimization.  Column  2 contains  the  parameter  ITYPE,  which  is  1 for  a 
complete  nozzle-base-boattail  assembly,  2 for  a nozzle,  or  3 for  a 
boattail.  The  remainder  of  the  card  (78  columns)  may  contain  any  alpha- 
numeric job  identification  data. 

The  following  cards  specify  the  remaining  input  data  in  free  format 
form,  read  in  through  NAMELIST  INFO.  Not  all  variables  need  to  be  de- 
fined for  each  run,  as  default  values  are  specified  in  subroutine  INPUTT 
for  most  of  the  variables. 

Each  sample  case  considers  the  same  thermodynamic  model,  which  is 
specified  by  the  default  values  in  the  program.  The  gas  in  the  nozzle 
is  a thermally  and  calorically  perfect  gas  with  y = 1.2  and  R = 60.0 
(ft-lbf)/(lbm-R).  Thus,  NGASN  = 1,  GAMN(l)  = 1.2,  and  RGASN(l)  = 60.0. 
The  gas  flowing  around  the  boattail  is  also  a thermally  and  calorically 
perfect  gas  with  y = 1.4  and  R = 53.3  (ft-lbf)/(lbm-R) . Thus,  NGASB  = 1, 
GAMB(l)  = 1.4,  and  RGASB(l)  = 53.3.  The  nozzle  gas  properties  are 
typical  of  high  temperature  exhaust  jet  gases,  and  the  boattail  gas 
properties  are  those  of  air.  The  stagnation  properties  are  specified: 
in  the  nozzle,  T0N  = 6000.0  R and  P0N  = 1000  psia,  and  for  the  boattail 
T0B  = 600.0  R and  P0B  = 12.0  psia. 

In  each  sample  case,  the  geometry  to  be  analyzed  is  discussed,  and 
the  running  time  on  the  CDC  6500  computer  is  given. 


1 . Sample  Case  1 

The  first  sample  case  illustrates  a nozzle  flow  field  analysis.  It 
considers  a nozzle  exhausting  into  a static  region  having  a pressure  of 
5.0  psia  (PAMB  = 5.0).  Thus,  in  column  1 of  the  title  card  I0PT  = 1 and 
in  column  2 of  the  title  card  ITYPE  = 2.  The  job  identification  state- 
ment "SAMPLE  CASE  1"  is  contained  in  the  remainder  of  the  title  card. 


A Kliegel  type  initial-value  line  is  specified  by  default  (IVLN  = 0). 
The  nozzle  geometry  is  specified  with  the  nozzle  upstream  radius  of 
curvature  RUN  - 3.0  in.,  the  nozzle  downstream  radius  of  curvature 
RDN  = 0.5  in.,  the  nozzle  tiiroat  radius  YTN  = 1.0  in.,  the  throat  attach- 
ment angle  THETTN  = 30.0  deg,  the  nozzle  exit  radius  YEN  = 5.0  in.,  and 
the  nozzle  length  XEN  = 10.0  in.  The  nozzle  wall  is  a second-order 
polynomial  contour  (IWALLN  =2).  All  of  the  parameters  are  the  program 
default  values. 

The  accuracy  of  the  analysis  depends  on  the  number  and  spacing  of 
the  initial-value  line  points.  In  this  sample  case,  11  points  are  speci- 
fied spanning  the  nozzle  throat  radius  (NPN  " 11). 

The  best  accuracy  is  obtained  when  these  points  are  spaced  geome- 
trically across  the  throat  with  the  ratio  of  the  distance  between  the 
two  points  at  the  centerline  to  the  distance  between  the  two  points  at  the 
wall  having  a value  larger  than  unity.  For  this  sample  case,  RATI0I  = 
5.0.  That  spacing  makes  the  change  in  property  values  between  points  of 
the  computational  mesh  roughly  comparable,  because  property  gradients  are 
greatest  near  the  wall.  A similar  procedure  is  employed  for  points  along 
the  nozzle  wall  on  the  initial  expansion  contour.  A wall  point  spacing 
compatible  with  the  initial-value  line  spacing  is  obtained  with  NPW  = 22 
and  RATI0W  = 5.0.  It  should  be  noted  that  the  number  of  points  in  the 
flow  field  solution  mesh,  and  hence  the  computation  time,  is  propor- 
tional to  the  square  of  NPN.  Large  amounts  of  computation  time  may  be 
saved  by  reducing  NPN,  although  a corresponding  decrease  in  accuracy  re- 
sults. The  above  parameters  are  all  specified  by  the  program  default 
values. 

The  parameters  T0LN  and  ICN  are  the  fractional  convergence  toler- 
ance and  the  maximum  number  of  applications  of  the  corrector,  respec- 
tively, in  the  method  of  characteristics  unit  processes.  Those  para- 
meters are  specified  by  their  default  values  (T0LN  = 0.001  and  ICN  = 2). 

Output  is  kept  brief  by  setting  KWRITN  = 1.  This  option  prints 
only  ini ti al -val ue  line  points,  axis  points,  wall  points,  and  the  one- 
dimensional approximation  to  the  solution  at  the  wall  points. 

The  following  list  presents  the  data  deck  for  sample  case  1.  Note 
that  only  KWRITN  and  PAMB  are  specified,  because  default  values  are  em- 
ployed for  all  of  the  other  parameters. 


12  SAMPLE  CASE  1 NOZZLE  ANALYSIS 
SINFO  KWRITN  = 1,  PAMB  = 5.0  $ 

A complete  listing  of  the  output  from  sample  case  1 is  presented  in 
Fig.  F-1  . The  first  page  contains  the  program  abstract  and  the  job 
title.  The  nozzle  flow  field  analysis  data  and  i ni ti al -val ue  line  are 
listed  on  the  second  page.  Column  headings  for  the  initial-value  line 
data  are  as  fol lows . 


69 


downstream  distance  from  the  nozzle  throat,  in. 
radial  distance  from  the  nozzle  axis,  in. 

Mach  number 
pressure,  psia 
velocity  magnitude,  ft/sec 
flow  angle,  deg. 
x-component  of  velocity,  ft/sec 
y-component  of  velocity,  ft/sec 
density,  Ibm/ft^ 
temperature,  R 

mass  flow  rate  between  the  axis  of  symmetry  and  the 
given  point,  Ibm/sec 

thrust  between  the  axis  of  symmetry  and  the  given  point, 

Ibf 

specific  impulse  for  the  flow  between  the  axis  of  synmetry 
and  the  given  point,  Ibf-sec/lbm 

The  initial-value  line  performance  parameters  are  also  listed. 

The  next  three  pages  of  output  presented  in  Fig.  F-1  contain  the 
results  of  the  method  of  characteristics  flow  field  analysis.  The  first 
of  these  three  pages  presents  the  wall  contour  specification  and  the 
column  headings  for  the  method  of  characteri sties  calculations.  J and  I 
are  the  characteristic  coordinates,  and  the  next  ten  columns  are  the 
same  as  described  above  for  the  initial-value  line  output.  Columns  13 
and  14,  headed  by  MDT/LRC  and  F/LRC,  respectively,  refer  to  the  inte- 
grated values  of  the  mass  flow  rate  and  the  thrust  along  the  left- 
running Mach  line  from  the  axis  of  symmetry  to  the  given  point.  These 
two  columns  represent  checks  on  the  accuracy  of  the  program.  The  mass 
flow  rate  at  any  wall  point  should  be  close  to  the  mass  flow  rate  at 
the  throat.  The  thrust  obtained  along  the  wall,  THRUST,  should  be  close 
to  the  value  obtained  along  the  left-running  characteristic,  F/LRC. 

The  specific  impulse,  ISP,  is  printed  in  the  next  column,  and  the  number 
of  applications  of  the  corrector,  IC,  or  the  specific  impulse  efficiency, 
EFF,  is  printed  in  the  last  column.  IC  is  printed  when  the  point  is  a 
flow  field  point  computed  by  any  of  the  method  of  characteristics  unit 
process,  and  EFF  is  printed  when  the  point  is  the  one-dimensional  approx- 
imation to  the  wall  point  solution.  When  0 is  printed  in  the  last 
column,  the  point  is  on  the  i ni tial -val ue  line.  A -2  in  the  last 
column  indicates  that  the  point  is  the  nozzle  exit  lip. 

Wall  points  are  located  at  the  end  of  every  left-running  Mach  line. 
However,  in  the  throat  region  of  the  nozzle,  the  inverse  wall  point 
method  is  employed.  That  procedure  arbitrarily  specifies  the  wall  point 
locations,  thus  controlling  the  spacing  and  enhancing  the  accuracy  of 
the  flow  field  analysis.  In  this  region  there  may  be  several  wall 
points  along  the  same  left-running  characteristic.  After  a wall  point 
is  printed,  the  one-dimensional  flow  approximation  to  the  solution  at 
the  wall  point  is  computed  and  printed  for  reference. 


X 

Y 
M 
P 
Q 
TH 
U 

V 
R 
T 

WD0T 

THRUST 

ISP 
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The  last  line  of  the  output  gives  the  total  number  of  characteris- 
tics that  crossed  inside  of  the  nozzle.  The  formation  of  oblique  shock 
waves  is  indicated  when  this  is  nonzero. 

Sample  case  1 required  5 seconds  of  central  processor  time  on  the 
CDC  6500  computer. 
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Figure  F-1.  Output  for  Sample  Case  1 (Continued). 


Figure  F-1.  Output  for  Sample  Case  1 (Continued). 
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Sample  case  2 illustrates  a boattail  flow  field  analysis.  Thus,  the 
title  cards  contains  13  in  tlie  first  two  columns.  A uniform  flow 
initial-value  line  is  specified  by  default  (IVLB  = 0).  The  wall  geometry 
is  specified  with  outer  radius  YTB  = 8.0  in.,  initial  expansion  contour 
radius  ROB  = 5.0  in.,  attachment  angle  THETTB  = 2.0  deg,  exit  radius 
YEB  = 7.0  in.,  and  length  XEB  = 10.0  in.  The  boattail  wall  is  the  second- 
order  polynomial  ( IWALLB  = 2).  All  of  these  parameters  are  the  default 
values  given  in  the  program.  For  this  sample  case,  KWRITB  = 2 for  abbre- 
viated output.  This  option  dictates  that  only  wall  points  are  printed 
out.  The  thrust  values  printed  in  the  last  column  are  obtained  by  inte- 
grating the  pressure  forces  acting  on  the  boattail  contour.  The  following 
list  presents  the  data  deck  for  sample  case  2. 

13  SAMPLE  CASE  2 BOATTAIL  ANALYSIS 
SINFO  KWRITB=2  $ 

Figure  F-2  presents  the  output  for  sample  case  2.  The  first  page 
presents  the  program  abstract  and  the  job  title.  The  second  page  pre- 
sents boattail  flow  field  analysis  data  and  the  i ni ti al -val ue  line.  The 
boattail  geometric  specification  is  given  at  the  top  of  the  third  page. 

The  third  and  fourth  pages  present  the  boattail  flow  field  solution. 

Sample  case  2 required  4 seconds  of  central  processor  time  on  the 
CDC  6500  computer. 
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Figure  F-2.  Output  for  Sample  Case 
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3.  Sample  Case  3 


Sample  case  3 illustrates  a complete  nozzle-base-boattai 1 flow 
field  analysis  for  the  combined  nozzle  and  boattail  contours  specified 
in  sample  cases  1 and  2.  Consequently,  11  is  placed  in  the  first  two 
columns  of  the  title  card.  The  total  thrust  on  the  assembly  is  the 
sum  of  the  thrusts  acting  on  the  nozzle,  the  boattail , and  the  base 
region.  For  this  sample  case,  KWRITN  = 2 and  KWRITB  = 2.  This  print 
option  dictates  that  wall  points  only  will  be  written  out.  The  fol- 
lowing list  presents  the  data  deck  for  sample  case  3. 


11  SAMPLE  CASE  3 NOZZLE-BASE-BOA TTAIL  ANALYSIS 

$INF0  KWRITN=2,  KWRITB=2,  YEB=5.2  $ 

Figure  F-3  presents  the  output  for  sample  case  3.  The  first  page 
of  the  output  lists  the  program  abstract  and  the  job  title  for  the  run. 
The  second  page  presents  the  nozzle  flow  field  analysis  data  and  initial- 
value  line  The  third  page  presents  the  nozzle  wall  geometric  specifi- 
cation, and  the  third  and  fourth  pages  present  the  nozzle  flow  field 
solution.  At  the  end  of  the  flow  field  data,  -2  appears  in  the  last 
column  denoting  the  nozzle  exit  lip  point.  The  number  of  internally 
crossing  characteristics  is  zero. 

The  three  pages  following  the  nozzle  flow  field  results  present 
the  boattail  flow  field  results.  The  first  page  lists  input  variables, 
initial -val ue  line  data,  and  geometric  parameters.  The  next  page 
presents  the  boattail  flow  field  solution.  For  this  sample  case, 

KWRITB  = 2,  thus,  only  wall  point  solutions  are  printed  out.  Each 
right-running  characteristic  begins  on  the  initial-value  line  and 
extends  to  the  boattail  wall,  where  the  thrust  is  computed.  This 
procedure  continues  until  a right-running  characteristic  passes  the 
exit  lip  point.  Then  the  inverse  wall  point  routine  computes  the 
boattail  exit  lip  point  properties. 

The  base  pressure  and  the  contribution  to  the  thrust  from  that 
pressure  are  then  calculated.  The  exit  lip  point  properties  of  the 
nozzle  and  the  boattail,  as  well  as  the  base  region  properties,  are 
printed  out  on  the  last  page.  The  total  thrust,  which  is  written  out 
on  the  last  line,  is  the  sum  of  the  nozzle,  boattail,  and  base  thrusts. 

Sample  case  3 required  13  seconds  of  central  processor  time  on 
the  CDC  6500  computer. 
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Figure  F-3.  Output  for  Sample  Case 
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Figure  F-3.  Output  for  Sample  Case  3 (Continued). 


Figure  F-3.  Output  for  Sample  Case  3 (Continued). 
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Figure  F-3.  Output  for  Sample  Case  3 (Continued) 


Figure  F-3.  Output  for  Sample  Case  3 (Concluded). 


4 . Sample  Case  4 


Sample  case  4 illustrates  a nozzle  optimization  using  Newton's 
method.  The  nozzle  gas  dynamic  model  and  geometry  are  specified  by 
the  program  default  valves  (see  Sample  Case  1).  To  reduce  the  output, 
KWRITN  = 3.  The  locations  of  the  initial-value  points  are  specified  by 
NPN  = 9 and  RATI0I  = 4.0.  The  locations  of  the  points  on  the  initial 
expansion  contour  are  specified  by  NPW  = 12  and  RATI0W  = 3.0.  The  ini- 
tial estimate  for  the  throat  attachment  angle  is  ANSTRT(l)  = 25.0.  To 
minimize  the  flow  field  output  during  the  optimization  procedure, 

KWRIT0  = 3.  The  nozzle  attachment  angle  perturbation  DA  = 2.0.  For 
a nozzle  optimization,  I0P  = 1 specifies  the  geometric  model  and 
N = 2 specifies  a two-dimensional  optimization.  The  following  list 
presents  the  data  deck  for  sample  case  4. 


22  SAMPLE  CASE  4 NOZZLE  OPTIMIZATION 
SINFO  KWRITN=3,  NPN=9,  RATI0I=4.O,  NPW=12,  RATI0W=3.O, 

ANSTRT(1)=25.0,  KWRIT0=3,  DA=2.0  $ 

The  first  page  of  the  output  prints  the  program  abstract  and  the 
job  title.  The  second  page  presents  the  input  data,  the  initial-value 
line,  and  the  initial-value  line  performance  parameters.  For  this  sample 
case,  KWRITN  = 3 and  KWRIT0  = 3 to  minimize  output.  A secondary-start 
line  is  computed  at  the  minimum  attachment  angle,  ANMIN(l)  = 20.0,  which 
is  printed  out  on  the  next  page  (see  the  discussion  in  Section  II).  The 
optimization  parameters  follow  on  the  next  page,  which  lists  the  type 
of  geometry  (I0P  =1),  print  option  (KWRIT0  = 3),  convergence  tolerance 
(T0L0  = 0.001),  optimization  variable  space  (N  = 2),  optimization  method 
(IMETH  = 3),  and  neximum  number  of  base  point  neves  (ITER0  = 10).  The 
initial  estimate  for  the  throat  attachment  angle  is  25.0  deg,  and  for 
the  exit  lip  radius  is  5.0  in.  These  values,  and  the  minimum  and  maximum 
values  of  those  variables,  are  listed  for  reference. 

The  next  three  pages  present  the  results  of  the  optimization.  On 
the  first  of  those  three  pages,  a heading  is  printed  out  specifying  the 
optimization  method.  The  next  three  lines  present  the  results  of  the 
nozzle  flow  field  analysis  for  the  initial  nozzle  contour.  The  output 
for  each  nozzle  flow  field  analysis  is  reduced  to  three  lines.  The 
first  line  specifies  the  nozzle  geometry,  the  second  line  presents  the 
flow  properties  at  the  nozzle  exit  lip  point,  and  the  third  line  is  the 
crossing  characteristics  count.  After  the  analysis  of  the  initial  nozzle 
contour,  an  optimization  base  point  data  line  is  written  out  stating  that 
for  the  base  point  (0  step),  THETA  = 25.0  deg,  YE  = 5.000  in.,  and  the 
nozzle  thrust  is  5644.65  Ibf. 

Newton's  method  requires  the  values  of  the  thrust  for  the  base  point 
and  five  neighboring  points.  The  nozzles  corresponding  to  these  five 
points  are  analyzed,  and  the  results  presented  in  the  next  five  groups 
of  three  lines  each.  From  that  data,  Newton's  method  constructs  the 
base  point  move  and  analyzes  the  corresponding  new  base  point  nozzle. 
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The  next  three  lines  present  the  results  of  that  analysis.  An  optimiza- 
tion base  point  data  line  is  then  written  out  presenting  the  base  point 
data  for  that  step  (1  step). 

The  above  procedure  is  repeated  for  a total  of  three  base  point 
moves,  at  which  time  the  optimization  procedure  has  converged.  The  last 
line  of  output  shows  that  the  optimum  nozzle  contour  has  a throat  attach- 
ment angle  of  32.729  deg,  an  exit  lip  point  radius  of  5.244  in.,  and  a 
thrust  of  5679.81  Ibf. 

Sample  case  4 required  80  seconds  of  central  processor  time  on 
the  CDC  6500  computer. 
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Figure  F-4.  Output  for  Sample  Case 


Figure  F-4.  Output  for  Sample  Case  4 (Continued). 


Figure  F-4.  Output  for  Sample  Case  4 (Continued). 


5 . S^cynp  1 e_  Case  5 

Sample  case  5 illustrates  a nozzle-base-boattai 1 assembly  optimi- 
zation. The  maximum  thrust  contour  is  approached  in  two  steps.  First, 
a two-dimensional  optimization  is  performed.  This  is  accomplished  by 
requiring  the  width  of  the  base  region  to  be  the  minimum  allowable  value 
(■Veb  " Yen  "yb)’  which  leaves  two  parameters  (0^^  yg^,)  uncon- 

strained. This  optimization  converges  quickly  because  of  the  relatively 
few  function  evaluations  required  for  a base  point  move.  After  conver- 
gence, the  optimization  problem  is  redefined  so  that  the  boattail  and 
nozzle  exit  lip  points  are  independent,  thus  yielding  a three-dimensional 
optimization  problem.  From  the  preliminary  two-dimensional  optimization, 
a good  initial  approximation  to  the  maximum  thrust  nozzle-base-boattail 
contour  is  available,  and  the  three-dimensional  optimization  converges  in 
one  pass. 

The  general  operating  conditions  considered  in  sample  case  5 are 
the  same  as  those  in  sample  case  3.  The  boattail  outer  radius  YTB  = 5.3 
in.,  the  radius  of  curvature  of  the  boattail  initial  expansion  contour 
is  the  program  default  value  RDB  = 5.0  in.,  and  the  boattail  contour 
is  conical  (IWALLB  = 1).  Output  is  minimized  by  setting  KWRITN,  KWRITB, 
and  KWRIT0  equal  to  3.  Computation  time  for  each  analysis  is  kept  small 
by  setting  NPN  = 9 and  NPB  = 40.  The  nozzle  has  12  inverse  wall  points 
(NPW  = 12)  spanning  the  minimum  attachment  angle  [ANMIN(l)  = 20.0  deg.]. 
The  initial -value  line  points  and  the  wall  points  are  spaced  according 
to  RATI0I  = 3.0,  RATI0W  = 3.0,  and  SPACE  = 0.6. 

Now  the  optimization  employs  Newton's  method  (IMETH  = 3)  to  calcu- 
late the  maximum  thrust  contours  to  within  a relative  tolerance  of 
T0L0  = 0.0001  or  within  a given  number  of  iterations  ITER0  = 10.  An 
initial  estimate  for  the  nozzle  geometry  is  given  by  specifying 
ANSTRT(l)  = 30.0  deg  and  ANSTRT(2)  = 4.0  in.  The  maximum  values  of 
6an>  yen>  ^ab’  Yeb  specified  as  ANMAX(l)  = 40.0  deg,  5.0  in., 

20.0  deg,  and  5.0  in.,  respectively.  The  minimum  values  of  those 
parameters  are  specified  as  the  program  default  values.  For  a nozzle- 
base-boattail  optimization,  I0P  = 3 specifies  the  type  of  geometry,  and 
N = 3 specifies  that  a three-dimensional  optimization  is  to  be  performed 
(i.e.,  the  unconstrained  geometric  variables  are  the  nozzle  throat  attach- 
ment angl  e Qan’  the  nozzle  exit  lip  radius  ygp,  and  the  boattail  exit 
lip  radius  yeb-  The  Gan  perturbation  DA  = 2.0  deg,  the  yen  perturbation 
DB  = 0.1  ytn>  ^nd  the  yg^  perturbation  DD  = 0.2  y^;f,.  The  minimum  annular 
base  width  DBAS  = 0,1  in. 

The  following  list  presents  the  data  deck  for  sample  case  5. 


21  SAMPLE  CASE  5 NOZZLE-BASE-BOATTAIL  OPTIMIZATION 
$INF0  KWRITN=3,  KWRITB=3,  KWRIT0=3,  YTB=5.3,  IWALLB=1, 
NPN=9,  NPW=12,  NPB=40,  RATI0I=3.O,  RATI0W=3.O,  SPACE=0.6, 
ANSTRT(1)=30.0,  4.0,  ANMAX(1 )=40.0,  5.0,  20.0,  5.0, 
IMETH=3,  ITER0=1O,  T0L0=O.OOO1,  I0P=3,  N=3, 

DA=2.0,  DB=0.1,  DBAS=0.1  $ 
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Figure  F-5  presents  the  output  for  sample  case  5.  The  first  page 
presents  the  program  abstract  and  the  job  title.  The  second  page  pre- 
sents the  nozzle  flow  field  analysis  data,  the  nozzle  initial-value 
line,  and  the  nozzle  i ni tial - val ue  line  perfoniiance  parameters.  The 
third  page  presents  the  nozzle  secondary  start  line.  The  optimization 
parameters,  the  initial  values  of  the  nozzle  geometric  parameters,  and 
the  minimum  and  maximum  values  of  the  nozzle  and  boattail  geometric 
parameters  are  presented  on  the  fourth  page.  The  fifth  page  presents 
the  boattail  flow  field  analysis  data  and  the  boattail  initial-value 
line. 


The  results  of  the  optimization  procedure  are  presented  on  pages 
six  to  ten  of  Fig.  F-5.  The  first  line  on  page  six  identifies  the 
optimization  procedure  as  Newton's  method.  The  next  group  of  seven 
lines  summarizes  the  results  of  the  combined  nozzle-base-boattail  flow 
field  analysis  for  the  initial  contour.  The  first  line  specifies  the 
nozzle  geometry,  the  second  line  present  the  flow  properties  at  the 
nozzle  exit  lip  point,  and  the  third  line  presents  the  crossing  character- 
istics counter.  The  fourth  line  presents  the  boattail  geometry,  the 
fifth  line  presents  the  flow  field  properties  at  the  boattail  exit 
lip  point,  and  the  six  line  present  the  base  region  pressure  and  thrust. 
Line  seven  presents  the  total  nozzle-base-boattail  thrust.  Each  flow 
field  analysis  for  a particular  nozzle-base-boattail  configuration 
generates  a corresponding  set  of  seven  lines  of  output. 

After  the  analysis  of  the  initial  nozzle-base-boattail  configuration, 
an  optimization  base  point  data  line  is  written  out  summarizing  the 
geometry  and  the  corresponding  thrust.  The  initial  contour  is  desig- 
nated the  0 step. 

The  two-dimensional  optimization,  in  which  the  boattail  exit  lip 
point  is  constrained  to  be  a fixed  distance  from  the  nozzle  exit  lip 
point,  then  commences.  This  procedure  is  analogous  to  the  nozzle  opti- 
mization discussed  in  sample  case  4.  Newton's  method  requires  the  values 
of  thrust  at  the  base  point  and  five  neighboring  points.  Those  five 
points  are  the  first  five  points  after  the  0 step  data  line.  From  that 
data,  Newton's  method  constructs  the  base  point  nwve  and  analyzes  the 
new  base  point  configuration.  The  sixth  set  of  data  after  the  0 step 
presents  the  results  of  that  analysis.  An  optimization  base  point  data 
line  is  then  written  out  presenting  the  base  point  data  for  that  step 
(1  step). 

The  above  constrained  boattail  optimization  is  repeated  for  a total 
of  three  base  point  moves.  The  final  results  are  printed  out  on  the 
eight  page  of  Fig.  F-7,  as  the  seven  lines  immediately  preceding  the 
optimization  data  line  for  step  3.  At  this  time,  a full  three-dimension- 
al optimization  is  initiated  in  which  the  conical  boattail  moves  inde- 
pendently of  the  nozzle  exit  lip  point.  The  initial  contour  for  the 
three-dimensional  optimization  is  chosen  to  be  the  nozzle  contour  ob- 
tained in  the  constrained  boattail  optimization  described  above, and  a 
boattail  exit  lip  point  moved  outward  from  the  nozzle  exit  lip  point  by 
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twice  the  minimum  separation  distance  Ay^.  That  contour  is  then 
analyzed  and  printed  out  as  the  next  seven  lines  of  output.  The  corre- 
sponding optimization  base  point  data  line  is  identified  as  step  0 for 
the  three-dimensional  optimization. 

For  the  three-dimensional  optimization,  Newton's  method  requires 
the  values  of  the  thrust  for  the  base  point  configuration  and  nine 
neighboring  points.  Those  results  are  presented  on  the  last  three 
pages  of  Fig.  F-5.  From  that  data,  Newton's  method  constructs  the  base 
point  move  and  analyzes  the  corresponding  new  base  point  configuration. 
The  final  group  of  seven  lines  on  the  last  page  of  Fig.  F-5  presents 
the  results  of  that  analysis.  The  last  line  on  Fig.  F-5  presents  the 
optimization  base  point  data  line  for  that  step,  step  1.  For  the  present 
sample  case,  the  results  of  step  1 satisfy  the  convergence  criteria  and 
the  program  stops . 

Sample  case  5 required  240  seconds  of  central  processor  time  on 
the  CDC  6500  computer. 
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Figure  F-5.  Output  for  Sample  Case 


Figure  F-5.  Output  for  Sample  Case  5 (Continued). 


Figure  F-5.  Output  for  Sample  Case  5 (Continued). 
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Figure  F-5.  Output  for  Sample  Case  5 (Continued). 
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